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ABSTRACT 


-This  thesis  lays  the  foundation  for  the  dynamic  stall 
Investigation  being  conducted  at  the  Fluid  Mechanics  Lab¬ 
oratory  at  NASA-Ames  Research  Center.  Using  existing  optical 
and  electrical  equipment,  a  new  dedicated  Mlcro-VAX  computer 
and  Labstar  software,  an  Indraft  transonic  wind  tunnel  and 
able  technicians  to  make  the  proper  Interface  hardware,  the 
project  came  together  In  a  new  test  facility  at  the  Fluid 
Mechanics  Laboratory.  The  goal  of  the  thesis  was  to  obtain 
both  qualitative  and  quantitative  Information  about  the  wake 
profiles  of  an  airfoil  in  steady  state  operations  at  varying 
angles  of  attack  and  tunnel  conditions.  To  accomplish  this 
task,  schlleren  photography  was  used  to  obtain  a  qualitative 
picture  of  the  flow  field.  With  this  information,  a  two  com¬ 
ponent  Laser  Doppler  Velocimeter  was  set  up  to  accurately 
measure  the  velocity  profiles  that  correspond  to  the  schlle¬ 
ren  photographs.  Once  this  preliminary  work  is  completed, 
the  same  apparatus  will  be  used  to  further  Investigate  the 
unsteady  dynamic  stall  phenomenon. 
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A.  DESCRIPTION  OF  THE  OVERALL  DYNAMIC  STALL  PROJECT 


Static  stall  Is  a  well  understood  and  documented  phenom¬ 
enon.  A  wing  or  rotor  blade  that  is  operating  at  a  high  <*, 
angle  of  attack,  will  experience  a  point  where  separated 
flow  will  cause  an  accelerated  decrease  In  lift  and  a  large 
rapid  increase  in  drag.  This  «*  at  which  the  wing  will  pro¬ 
duce  the  maximum  lift,  C]rnax,  Is  a  very  important  aspect 
of  airfoil  performance.  This  stall  a  determines  the  stall 
speed  of  the  aircraft.  Considerable  effort  has  been  expended 
to  find  ways  to  decrease  the  stall  speed  by  Increasing  the 
C]max-  The  majority  of  this  effort  has  been  directed  to¬ 
ward  a  gradual  approach  to  the  stall  «,  in  other  words, 
static  stall  research. 

The  wing  can  be  made  to  continue  to  produce  lift  past 
the  static  stall  a  if  the  wing  angle  of  attack  is  In¬ 
creased  rapidly  past  the  stall  <*.  The  phenomenon  of 
dynamic  stall  is  manifested  in  the  rapid  decrease  of  the 
lift  at  the  higher  stall  a.  There  are  four  basic  stages  of 
dynamic  stall  that  have  been  documented  [Ref.  1): 

1.  Trai 1 1 ng-edge  flow  reversal,  progressing  forward  along 
the  airfoi 1 . 

2.  Formation,  growth  and  shedding  of  the  vortex. 

3.  Full  stal 1 . 

4.  Boundary  layer  reattachment. 

The  understanding  of  this  phenomenon  will  be  beneficial  to 


the  helicopter,  where  constantly  oscillating  rotor  blades 
are  consistently  In  the  dynamic  stall  region  throughout 
their  revolution,  and  jet  fighters,  where  agile  and  rapid 
maneuvering  through  the  stall  a  is  required  for  surviva¬ 
bility.  By  quantifying  the  velocity  field  around  and  behind 
the  airfoil,  a  better  understanding  of  the  dynamics  of  this 
phenomenon  will  be  gained.  This  new  analytical  information 
will  provide  valuable  Information  to  the  design  engineer  to 
expand  the  flight  envelope  of  future  fighter  aircraft  and 
improve  the  performance  of  the  helicopter  rotor  system. 

B.  THESIS  GOALS 

This  thesis  lays  the  foundation  for  the  dynamic  stall 
investigation  project  being  conducted  at  the  Fluid  Mechanics 
Laboratory  at  NASA-Ames  Research  Center.  Using  existing 
optical  and  electrical  equipment,  a  new  dedicated  Micro-VAX 
computer,  new  Labstar  software,  an  existing  transonic  in¬ 
draft  wind  tunnel  and  able  technicians  to  make  the  proper 
interfaces,  the  project  came  together  in  a  new  test  cell  of 
the  Fluid  Mechanics  Laboratory.  The  goal  of  this  thesis  was 
to  first  obtain  stroboscopic  schlieren  photographs  of  a  sta¬ 
tionary  airfoil  at  different  a's  to  obtain  qualitative 
data  on  the  flow  past  the  airfoil  and  the  size,  shape  and 
extent  of  the  wake  and  then  set  up  a  Laser  Doppler  Velocl- 


meter  to  obtain  quantitative  velocity  data  in  order  to  con¬ 
struct  an  accurate  velocity  profile  of  the  flow  field. 


C.  SCHLIEREN  AND  LASER  DOPPLER  VELOCIMETRY  BASICS 
1 .  Sch 1 leren 

Many  optica!  flow  visualization  techniques  have  been 
developed  for  gas  flows.  Each  of  these  techniques  are  sensi¬ 
tive  to  variations  In  density  and  derivatives  of  density  In 
the  flow.  Interferometry  Is  sensitive  to  the  density  of  the 
gas  flow;  f.  Schl leren  is  sensitive  to  the  density  grad¬ 
ient  or  the  change  In  density  with  respect  to  distance; 
if  sty .  Shadowgraphy  Is  sensitive  to  the  second  derivative 
of  the  density;  i'Zf'i  y2.  All  the  techniques  rely  on 
the  fact  that  the  speed  of  light  Is  a  function  of  the  local 
index  of  refraction  in  that  medium.  In  gases,  the  Index  of 
refraction  Is  a  function  of  density. 

The  schl leren,  German  for  streak,  technique  allows 
one  to  view,  qualitatively,  the  density  gradient  In  a  given 
flow  field.  This  visualization  system  Is  considerably  sim¬ 
pler  and  cheaper  to  use  than  other  techniques.  Therefore,  It 
was  chosen  as  the  preliminary  technique  to  obtain  qualita¬ 
tive  information  about  the  air  flow  around  and  In  the  wake 
of  the  airfoil.  The  schl leren  setup  Involves  using  a  line 
source  of  light  and  projecting  It  through  the  test  section 
where  the  flow  field  alters  the  path  of  the  light.  The  line 
source  Is  usually  a  monochromatic  gas  discharge  tube.  The 
light  Is  projected  through  a  narrow  slit  to  provide  a  line 
of  light  of  uniform  Illumination.  The  dimensions  of  the  slit 
control  the  definition  of  the  image  seen  on  the  screen.  A 
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high  quality  concave  mirror  Is  usually  used  to  direct  the 
light  from  the  source  through  the  test  section.  On  the  other 
side  of  the  test  section,  another  similar  mirror  Is  used  to 
reflect  the  light  to  a  viewing  screen  or  camera.  Each  point 
of  the  line  source  will  illuminate  each  point  in  the  test 
section  if  the  source  is  collimated  and  of  equal  optical 
distance  on  both  sides  of  the  test  section.  The  light  from 
the  second  mirror  is  focused  onto  a  knife  edge.  The  place¬ 
ment  and  adjustment  of  the  knife  edge  Is  very  critical  as 
the  sensitivity  of  the  entire  system  depends  on  It.  The 
knife  edge  Is  placed  In  the  plane  of  the  source  Image  and 
interferes  with  a  portion  of  this  Image.  A  portion  of  the 
light  from  each  point  in  the  test  section  is  blocked  which 
decreases  the  light  impinging  on  the  viewing  screen.  The 
density  gradients  encountered  In  the  test  section  normal  to 
the  knife  edge  cause  a  refraction  of  light  allowing  more 
light  from  that  portion  of  the  test  section  to  pass  over  the 
knife  edge  and  Illuminate  the  viewing  screen.  By  rotating 
both  the  line  source  and  the  knife  edge  both  directional 
derivatives  of  density  can  be  viewed.  Also,  through  the 
lateral  adjustment  of  this  knife  edge  at  the  source  image, 
a  fine  control  over  the  degree  of  detail  resulting  from  this 
technique  can  be  exercised. 

The  viewing  screen  can  be  any  Image  recording 
device.  Polaroid  cameras  can  be  set  up  for  Instantaneous 
results  or  35  mm  cameras  can  be  mounted  at  the  focal  length 
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of  the  lens  to  obtain  numerous  photographs  during  the 
testing.  Video  cameras  can  also  be  used. 

2.  Laser  Doppler  Veloclmetrv 

The  techniques  for  measuring  velocity  in  the  past 
used  sensors,  which  need  to  be  introduced  into  the  flow. 
This  necessarily  meant  Intruding  Into  and  affecting  the 
flow.  During  the  last  two  decades,  a  new  technique.  Laser 
Doppler  Veloclmetry  or  LDV,  has  been  available.  LDV's  are 
Indeed  a  better  technique.  It  Is  nonlntrusive  and  Is  capable 
of  directly  measuring  velocities  of  particles  in  the  flow. 
Laser  Doppler  Veloclmetry  has  been  very  successful  In  meas¬ 
uring  flow  velocities  in  turbulent  and  reverse  flows.  It  has 
been  successfully  used  for  the  past  two  decades  to  obtain 
velocity  measurements  accurately. 

Laser  Doppler  Veloclmetry  uses  two  coherent!  light 
beams  from  the  same  source,  which  Intersect  to  form  the 
measurement  volume.  The  two  light  beams  form  an  interference 
fringe  system  that  Is  crossed  by  particles  moving  with  the 
flow.  The  light  scattered  by  the  particles  Is  Doppler 
shlfted2.  The  frequency  of  this  scattered  light  is 
directly  proportional  to  the  particle  velocity.  With 


1.  Coherent  Light  =  light  in  which  the  electromagnetic 
waves  maintain  a  fixed  phase  relationship.  (Random  House 
College  Dictionary) 

2.  Doppler  Shift  =  the  magnitude  of  the  apparent  change 
in  frequency  resulting  from  the  relative  motion  of  the 
light  source  and  the  receiver.  (Random  House  College  Diet) 
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suitable  selection  of  flow  particles,  known  also  as  seeds, 
the  fluid  velocities  can  be  measured  from  the  particle 
velocities.  Using  photodetectors,  the  light  scattered  by  the 
particles  Is  converted  into  electrical  signals.  This  signal, 
after  processing,  provides  frequency  information  from  which 
the  velocity  can  be  determined.  The  measurement  of  two 
components  of  the  velocity  can  be  obtained  by  setting  up  two 
orthogonal  systems  with  laser  beams  of  two  colors.  In  the 
same  optics  train,  thus  identifying  each  component  with  a 
particular  color. 


A.  SCHLIEREN  SETUP 


The  schlleren  system  for  this  experiment  was  set  up 
using  both  a  continuous  light  source  and  a  stroboscopic 
source.  The  continuous  source  was  used  to  provide  a  time 
averaged  schlleren  photograph  and  also  to  allow  use  of  a 
video  tape  recorder  to  record  the  time  history  of  the  flow. 
To  provide  Instantaneous  schlleren  photographs  of  the  flow 
field,  the  stroboscope  was  used.  See  Figure  2.1  for  the 
schematic  of  the  setup. 

To  obtain  the  35  mm  photographs  of  the  flow  field,  a  US 
Scientific  Strobrlte  Stroboscope  was  used.  The  strobe  had 
five  settings,  20  Hz  to  1  KHz,  to  adjust  the  duration  of  the 
strobe.  The  70  Hz  setting  proved  to  be  the  optimal  setting. 
A  mechanically  adjustable  slit  In  front  of  the  strobe  was 
used  to  produce  a  line  source  of  light  needed  for  the  tech¬ 
nique  to  operate  effectively.  A  large  concave  mirror  with  a 
305  cm  focal  length  reflects  this  light  sheet  towards  the 
test  section.  The  test  section  had  optical  access  through 
both  sides  of  the  tunnel.  The  test  section's  dimensions  were 
13  cm  by  25  cm  and  was  constructed  out  of  2.5  cm  thick 
Plexiglas.  The  light  passing  through  the  test  section  was 
then  reflected  by  an  Identical  mirror  located  opposite  the 
first  mirror  on  the  other  side  of  the  test  section  as  shown 
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Figure  2.1  Schlleren  Schematic 

In  Figure  2.1.  Light  from  the  second  mirror  was  focussed  on 
a  vertical  knife  edge  whose  position  could  be  adjusted  to 
get  a  highly  sensitive  system.  The  knife  edge  was  oriented 
vertically  to  obtain  the  visual  signature  of  the  density 
gradient  In  the  vertical  plane  of  the  test  section.  After 
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the  knife  edge,  a  200  mm  lens  was. used  to  focus  the  image 
onto  the  photographic  device. 

Since  the  test  section  is  25  cm  wide,  the  photographs 
obtained  from  the  system  actually  represent  density  gradient 
effects  averaged  over  the  test  section  width,  as  Is  always 
the  case  with  the  schlieren  technique. 

A  Nikon  35  mm  camera  body  with  motordrive  (but  no  lens) 
was  used  to  record  the  schlieren  Images.  The  camera  was 
placed  at  the  focal  point  of  the  200  mm  lens.  Fine  adjust¬ 
ments  to  focus  the  image  in  the  camera  were  made  by  a  micro¬ 
meter  slide  on  the  camera  mounting  mechanism.  Pictures  were 
taken  by  opening  the  shutter  with  a  remote  shutter  release 
cable,  before  the  strobe  was  fired  remotely  and  then  closing 
the  shutter.  With  a  little  practice,  this  procedure  was 
quite  rapid.  Photographs  were  taken  at  two  different  Mach 
numbers,  0.30  and  0.45,  and  at  three  different  <*'s;  0°, 
7.7°,  15.5°.  See  Figures  A.l  through  A. 6,  Appendix  A, 
for  examples  of  schlieren  photographs.  This  series  of  photo¬ 
graphs  was  studied  to  analyze  the  flow  qualitatively  and  to 
determine  the  locations  in  the  wake  for  obtaining  the  veloc¬ 
ity  data  for  the  two  different  Mach  numbers. 

B.  LASER  DOPPLER  VELOCIMETER  SETUP 

This  section  will  discuss  the  specific  LDV  system  setup 
in  test  cell  number  four  at  the  Fluid  Mechanics  Laboratory. 
The  following  sections  will  be  discussed:  transmitting 


optics,  receiving  optics  and  photodetectors,  test  section, 
signal  processing  and  data  processing.  Figure  2.2  shows  a 
photograph  of  the  transmitting  optics  and  Figure  2.3  shows  a 
photograph  of  the  receiving  side.  Individual  schematics  will 
be  discussed  In  detail  within  the  following  sections. 

The  laser  used  Is  a  Spectra  Physics,  Model  164,  5  Watt 
Argon-Ion  Laser  capable  of  outputting  light  from  0.351  to 
0.528  microns  and  up  to  twelve  discrete  wavelengths  simul¬ 
taneously  IRef.  2:p.  231.  By  choosing  the  blue  beam,  0.4880 
microns,  and  the  green  beam,  0.5145  microns,  for  the  two 
component  Laser  Doppler  Veloclmetry  System,  nearly  75%  of 
the  Initial  output  power  could  be  used.  This  standard 
practice  was  adapted  in  this  case. 

1 .  Transmitting  Optics 

Owing  to  traversing  requirements,  the  laser  was 
higher  on  the  optics  table  than  the  10.8  cm  optical  height 
of  the  transmitting  optics.  A  set  of  two  mirrors  at  45° 
angles  to  the  direction  of  the  incident  light  were  used  to 
lower  the  beam  to  the  require  height.  Figure  2.4  shows  a  de¬ 
tailed  schematic  of  the  entire  transmitting  optics  system. 

The  vertically  polarized  beam  exiting  the  multiline 
laser,  after  being  directed  by  mirrors,  is  passed  through  a 
beam  collimator,  part  of  the  TSI  Model  9106  set.  Incident 
light  Is  collimated  If  the  beam  diameter  does  not  change  as 
It  passes  through  space.  In  reality,  however.  It  is  suffic¬ 
ient  if  the  collimator  can  insure  this  quality  through  the 
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Optical  Components 

1- Green  Polarization  Rotator 

2- Green  Beam  Splitter 

3- Green  Bragg  Cell 

4- Blue  Beam  Displacer 

5- Blue  Polarization  Rotator 

6- Blue  Beam  Splitter 

7- Blue  Bragg  Cell 

8- Green  Steering  Wedges 

9- Blue  Steering  Wedges 

10- Beam  Stops  (Both  Beams) 

11- Beam  Spacers  (Both  Beams) 


Green  Mirror 
Blue  Mirror  I 


Beam  Lowering  Mirrors 

Collimator 
Polarization  Rotator 

Half-Wave  attenuator' 

Mirror 

Dispersion  Prism/Color  Separator 


Figure  2.4  Transmitting  Optics  Schematic 


relatively  short  optical  path  of  the  system.  Proper  col  11- 
matlon  Is  said  to  be  achieved  when  the  waist1  of  the  beam 
coincides  with  the  focal  point  of  the  transmitting  lens.  By 
sliding  the  collimator  sleeve  In  and  out  and  observing  the 
beam  "spots"  on  a  distant  wall,  the  collimator  can  be  set. 


1.  Waist  =  the  point  on  a  beam  where  the  diameter  is 
the  sma I  1 est . 
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From  the  collimator,  the  beam  passes  through  a  po¬ 
larization  rotator  and  an  attenuator.  The  polarization  ro¬ 
tator  changes  the  polarization  of  the  Incoming  laser  beam  to 
conform  to  the  polarization  requirements  of  the  next  optical 
component,  a  dispersion  prism.  The  attenuator  works  in 
tandem  with  the  rotator  to  insure  passage  of  only  the  horiz¬ 
ontally  polarized  light  through  these  components.  The  dis¬ 
persion  prism  separates  the  wavelengths,  or  colors,  of  the 
light.  The  two  primary  colors  of  this  laser,  green  (0.5145 
microns)  and  blue  <0.4880  microns),  are  used  to  form  the  two 
component  system.  A  set  of  mirrors,  TSI  Model  9107,  Is  used 
to  direct  each  of  the  beams  to  the  appropriate  axis  for  the 
transmitting  optics  as  shown  in  Figure  2.4.  As  can  be  seen 
from  the  figure,  the  green  beam  Is  directed  down  the  center- 
line  of  the  optics  and  the  blue  is  directed  25  mm  to  the 
right  of  the  green. 

The  centerline,  green,  beam  passes  through  a  polar¬ 
ization  rotator,  TSI  Model  9102,  to  rotate  the  beam  polari¬ 
zation  90°,  then  a  beam  splitter,  TSI  Model  9115.  The 
polarization  change  Is  necessary  to  Insure  proper  splitting 
Into  two  beams  from  the  incident  light.  The  beam  splitter 
divides  the  Incident  beam  Into  two  beams  of  equal  Intensity 
and  parallel  to  the  optical  axis,  spaced  50  mm  apart. 


Since  the  flow  field  to  be  investigated  will  have  a 


turbulent  wake  with  possible  separated  flow  regions,  a  freq¬ 


uency  shift  system  is  used  for  each  of  the  two  colored 


beams.  Upshlftlng  of  the  frequency  of  the  laser  light  prior 


to  entering  the  test  section  will  allow  detection  of  the 


direction  the  particle  traverses  through  the  probe  volume. 


The  shift  results  in  a  moving  fringe  pattern  that  will  gen¬ 


erate  a  photodetector  signal  equal  to  the  shift  frequency 


when  measuring  a  particle  with  zero  velocity.  This  system 


allows  a  greater  range  of  measurements  of  negative  veloc¬ 


ities  and,  thus,  is  required  for  this  system  to  obtain  valid 


measurements. 


The  pair  of  green  beams  now  enter  a  frequency  shift¬ 


ing  device  called  the  Bragg  cell.  This  device  is  an  acousto¬ 


optic  unit  that  upshifts  the  frequency  of  an  incident  beam, 


The  shift  is  set  at  40  MHz.  This  device  is  placed  in  the 


optical  path  of  one  of  the  two  beams  forming  the  probe  vol¬ 


ume  for  that  color.  The  outgoing  beam  will  be  shifted  by  a 


40  MHz  shift  relative  to  the  other  beam.  The  Bragg  cell 


operates  by  deflecting  one  beam  in  a  fan-like  pattern  in 


the  plane  intersecting  both  incident  beams.  Each  separate 


deflected  beam  has  a  frequency  that  is  a  multiple  of  40  MHz 


from  the  incident  beam.  The  Bragg  cell  should  be  adjusted 


to  shift  the  light  intensity  to  the  first  order  beam.  This 


,s  accomplished  by  tilting  the  Bragg  cell  in  its  mount 


This  first  order  beam  should  have  an  intensity  greater  than 
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90%  of  the  Incident  light  beam  [Ref.  3:p.  1481.  An  optical 
wedge  Is  included  in  this  device  at  the  output  end  of  the 
Bragg  Cell  to  steer  the  deflected,  first  order  beam  parallel 
to  the  unshlfted  beam.  The  result  of  this  Bragg  shifting  Is 
an  interference  pattern  In  the  probe  volume  that  will  move 
relative  to  a  fixed  particle.  This  moving  fringe  pattern  can 
be  used  to  differentiate  the  direction  of  motion  through  the 
probe  volume.  Particles  moving  with  the  direction  of  motion 
of  the  fringes  gives  a  lower  frequency  and  against  the 
fringe  movement  gives  a  higher  frequency  as  shown  in  the 
following  equation.  The  particle  moving  through  the  probe 
volume  will  produce  a  signal  whose  frequency  Is 

^signal  *  'fghlft  +/“  fDoppler' 

where  the  sign  is  +  for  particles  moving  in  the  opposite 
direction  to  the  fringe  movement  [Ref.  2:p.  67]. 

For  optic  tube  support  purposes  a  TSI ,  Model  9176, 
Support  is  used.  This  device  also  serves  to  increase  the 
effectiveness  of  the  beam  stops  of  the  blue  beam  by  inc¬ 
reasing  the  deflected  beam  separation  distance. 

The  off-axis  blue  beam  Is  now  brought  to  the  center- 
line  using  a  beam  displacer,  TSI  Model  9174.  This  model  em¬ 
ploys  a  rhomboid  prism.  It  displaces  the  beam  a  fixed.  25 
mm,  distance  to  the  centerline  of  the  optics  while  keeping 
the  beam  parallel  to  the  axis.  The  blue  beam  then  passes 
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through  the  same  series  of  optical  components  as  the  green 
beams:  polarization  rotator,  beam  splitter  and  Bragg  cell. 

The  next  two  elements  In  the  optics  train  are  beam 
steering  modules,  TSI  Model  9175,  one  for  each  of  the 
colors.  These  modules  allow  for  steering  of  the  shifted  beam 
to  Insure  proper  beams  crossing. 

The  Bragg  cel!  produces  several  beams,  of  which 
only  the  first  harmonic,  the  beam  that  Is  shifted  by  40 
MHz,  Is  Important.  All  the  other  unwanted  beams  are  blocked 
off  by  means  of  a  beam  stop.  In  the  present  case,  a  beam 
stop  was  used  on  each  of  the  two  shifted  beams. 

Following  the  beam  stops,  beam  spacers,  TSI  Model 
9113-13,  are  used  to  for  each  set  of  beams  to  reduce  the 
distance  between  the  beam  pairs  forming  the  probe  volume. 
The  spacing  was  reduced  to  13  mm  by  using  these  prisms; 
while  maintaining  polarity,  coherence,  beam  diameter  and 
parallel  optical  paths.  The  primary  use  of  these  elements  Is 
to  control  the  fringe  spacing  In  the  probe  volume  by 
changing  the  angle  of  intersection  of  the  beam. 
Additionally,  however,  they  permit  transmission  of  the  beams 
Into  the  test  section  by  requiring  only  small  mirrors  and 
optic  components. 

The  four  beams  from  the  optic  tube  are  directed  to 
the  focusing  lens  with  five  mirrors.  The  first  two  move  the 
beams  from  10.8  cm  to  23  cm  above  the  optical  table  and  the 
next  two  reflect  the  beams  back  180°  toward  the  last 


mirror  In  this  series.  This  mirror  directs  the  four  beams 
vertically  toward  the  transmitting  optics  lens.  The  four 
beams  then  enter  the  transmitting  optics  lens  to  focus  the 
four  beams  to  a  point.  A  final  mirror  is  used  to  direct  this 
focused  beam  to  the  test  section.  The  focusing  lens  and  the 
last  three  mirrors  are  mounted  on  a  traversing  mechanism  to 
enable  traversing  of  the  focused  beams  in  any  of  the  three 
orthogonal  directions  as  shown  in  Figure  2.6. 

2.  Receiving  Optics  and  Photodetector 

All  the  receiving  optics  are  mounted  on  a  traversing 
mechanism  which  Is  slaved  to  the  traversing  mechanism  on  the 
transmitting  side.  The  two  sides  were  aligned  to  within  0.01 
mm  and  could  be  moved  together.  A  Taskmaster  controller  is 
used  to  move  the  two  traverses  both  manually  and  through  the 
use  of  computer  programs.  See  Figure  2.5. 

The  receiving  optics  are  used  to  capture  the  light 
scattered  by  the  particles  crossing  the  probe  volume  and 
focus  this  light  onto  photodetectors.  It  consists  of  several 
components;  that  are  mounted  on  a  hardware  arm  connected  to 
the  traversing  mechanism.  The  scattered  light  Is  collected 
by  a  series  of  three  lenses  that  collimate  the  light  and 
direct  it  toward  the  polarization  color  separator.  The  last 
lens  acts  as  a  spatial  filter  In  addition  to  being  a  beam 
collimator.  Its  function  Is  to  allow  light  focused  on  the 
lens  to  pass  through  while  blocking  all  the  rest  of  the 
light.  It  is  used  to  reduce  the  stray  light  scattered  from 
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arouad  the  probe  volume,  as  Is  the  case  In  measurements 
close  to  a  solid  surface. 

The  polarization  beam  separator  separates  the  two 
wavelengths  present  In  the  scattered  light  Into  Its  con¬ 
stituent  colors;  green  and  blue.  This  splitter  reflects 
light  with  horizontal  polarization  upward  while  transmitting 
the  vertically  polarized  light  through  it.  This  Insures  sep¬ 
aration  of  the  two  primary  colors  in  the  scattered  light. 

Two  notch  filters  are  used  at  this  stage,  on  each  of 
the  colors  to  eliminate  residual  color  contamination.  The 
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Figure  2.5  Receiving  Optics  Schematic 


separated  scattered  light  Is  now  focused  on  to  the  photo¬ 
detector  aperture  with  a  lens. 

The  photodetector  converts  the  light  signal  Into  an 
electric  signal.  The  photodetector  used  is  a  photomultiplier 
tube  or  PMT.  The  advantages  of  the  PMT  over  the  other  types 
are  low  noise,  high  gain  and  broad  bandwidth  while  Its 
fragility  and  the  need  for  a  high  voltage  source  are  its 
disadvantages  [Ref.  2:p.  451. 

3.  Test  Section 

The  transonic  test  section  is  formed  by  Inserting 
upper  and  lower  ventilated  walls  In  a  subsonic  Indraft  wind 
tunnel.  The  test  section's  dimensions  are  13  cm  by  25  cm. 


Figure  2.6  Transmitting  Side  of  Test  Section 


A  7.5  cm  chord  NACA  0012  airfoil  was  mounted  between  the 
side  walls  at  the  centerline.  This  airfoil  is  attached  at 
quarter  chord  to  a  external  device  to  oscillate  the  airfoil 
through  various  angles  of  attack.  The  side  walls  of  the  test 
section  are  made  of  2.5  cm  thick  plexiglas  to  allow  complete 
and  unobstructed  optical  access  to  the  airfoil.  See  Figure 
2.6.  The  mechanism  to  mount  the  model  partially  blocks  the 
views  of  the  leading  edge  of  the  model.  However,  this  was 
not  a  limitation  for  the  present  study. 

A.  Signal  Processing 

The  signal  from  the  photodetector  must  be  processed 
to  obtain  the  frequency  information  in  the  scattered  light. 
Figure  2.7  shows  the  block  diagram  of  the  signal  processing 
components.  The  signal  from  the  PMT  is  amplified  by  an  EIN 
Model  403LA,  linear  amplifier  to  boost  signal  strength  by  37 
dB  to  transmit  the  signal  over  a  long  distance  and  to  be  of 
appropriate  strength  for  signal  processing. 

A  downmixing  circuit  is  used  next  to  decrease  the 
effective  frequency  shift  from  40  MHz  to  around  10  Mhz  to 
improve  the  counter  resolution.  The  two  elements  of  the 
downmlxer  are  a  Tektronix  SG503  Levelled  Sine  Wave  Generator 
and  a  Hewlett-Packard,  Model  10534A,  Passive  Diode  Mixer. 
The  expected  Doppler  frequency  should  be  estimated  to  prop¬ 
erly  set  the  mixer  frequency.  Using  the  range  of  Mach  num¬ 
bers  expected  as  -0.2  to  +0.5,  the  free  stream  velocity 


range  can  be  calculated  using  the  set  of  equations  in 
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be  1.6  to  13.6  MHz.  The  mixing  frequency  should  be  adjusted 
whenever  the  tunnel  conditions  change  or  In  attempts  to 
obtain  better  counter  resolution. 

The  next  signal  processing  unit  encountered  is  the 
Series  3000  Macrodyne  LDV  Processor.  This  device  has  several 
functions:  variable  high  and  low  pass  filters.  Input  and 
output  ports,  threshold  and  overload  adjustments  and  freq¬ 
uency  counter  with  fringe  passing  selection  switch.  Both  the 
low  pass  filters,  to  remove  the  high  frequency  noise,  and 
the  high  pass  filters,  to  reduce  the  low  frequency  pedestal 
from  the  signal,  are  adjustable  to  provide  the  optimal  fil¬ 
tering  of  the  signal.  Once  the  signal  has  been  thus  condi¬ 
tioned,  a  frequency  counter  is  used  to  measure  the  signal 
frequency  and  convert  It  to  digital  time  data  that  can  be 
sent  to  the  data  processor.  The  counter  operates  on  two 
different  systems  to  validate  the  Doppler  signal:  a  three 
level  detection  scheme  and  a  5/8  or  10/16  comparison.  The 
three  level  detection  scheme  shown  In  Figure  2.0  must  be  met 
before  a  signal  is  validated.  A  threshold  level  Is  set  on 
the  counter  to  be  Just  above  the  level  of  the  steady  state 
noise.  The  signal  must  pass  the  following  gates  to  be  ac¬ 
cepted:  <1>  exceeds  the  threshold  setting  <2>  crosses  zero 
and  <3)  exceeds  a  negative  threshold.  To  pass  this  check, 
the  selected  number  of  fringe  crossings  must  sequence 
through  each  check  without  error.  This  detection  scheme  Is 
notably  effective  in  noise  discrimination.  Once  a  signal  is 
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received  and  passes  the  first  checK,  the  counter  starts  to 
measure  the  Doppler  period  by  opening  a  time  gate  and  stop¬ 
ping  at  the  switch  set  number  of  cycles.  This  setting  of 
cycles  is  able  to  be  selected  from  a  switch  on  the  counter 
control  panel.  Validated  data  from  the  switch  being  in  the 
10/16  cycles  position  tend  to  be  better  for  strong  Doppler 
bursts  because  the  greatest  number  of  cycles  measured  tends 
to  reduce  the  effects  of  noise  in  the  signal.  However,  this 
setting  also  restricts  the  motion  of  the  particle  through 
the  probe  volume  and  thus  Introduces  artificial  biases  of 
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Figure  2.8  Three  Level  Counter  Detection  Scheme 


the  data.  The  particle  must  traverse,  the  probe  volume  at 
Its  largest  length  to  Insure  validated  data  for  16  fringe 
crossings  For  this  reason,  present  experiments  were  carried 
out  with  the  switch  set  to  the  5/8  position. 

5.  Data  Processing  Hardware 

Once  the  data  Is  received  from  the  Macrodyne  signal 
processors.  It  must  be  multiplexed  and  transferred  to  the 
computer  for  velocity  calculations.  The  NASA  LDV  Computer 
Interface  Is  now  used  to  control  the  multiplexing  of  data, 
sequencing  of  data,  provide  a  digital  display  to  track 
successful  counts  of  the  Doppler  frequency  and,  generally, 
control  the  flow  of  data  between  the  signal  processor  and 
the  Mlcro-VAX  computer.  Since  a  two  component  Laser  Doppler 
Veloclmeter  Is  used  In  this  experiment,  the  multiplexer 
sends  four  channels  of  data  for  each  valid  Doppler  burst. 
The  first  data  word  contains  the  time  Information,  the 
second  contains  any  status  Information  from  the  signal  and 
the  last  two  words  contain  the  frequency  data  from  each 
component.  Due  to  the  limited  nature  of  this  experiment, 
the  first  two  words  of  data  are  not  useful  In  further 
velocity  calculations.  Therefore,  Just  the  frequencies  from 
each  particle  caused  Doppler  burst  are  used  to  calculate  the 
respective  particle  velocities. 

The  Mlcro-VAX  computer  stores  all  acquired  data  In  a 
buffer  set  up  by  acquisition  programs.  The  next  step  after 
data  acquisition  Is  to  calculate  the  velocities.  These 


computer  programs  reduce  the  massive  amounts  of  raw  data 
obtained  Into  meaningful  velocity  data  that  can  be  output  In 
tabular  form  and  graphical  form. 

The  following  is  a  simplified  path  followed  by  the 
numerous  computer  programs  to  reduce  the  data.  At  least  1000 
samples  from  each  component  are  required  to  obtain  a  statis¬ 
tically  steady  sample  frequency.  This  frequency  data  must 
then  be  converted  to  the  true  Doppler  frequency  by  the  fol¬ 
lowing  equation: 

fd  “  fs  “  fB  ♦  fM 

where  f^  Is  the  true  frequency  resulting  from  the  particle 
traversing  the  probe  volume*  fs  is  the  signal  frequency 
measured  by  the  counter,  fg  Is  the  Bragg  shift  frequency 
of  40  MHz  and  fM  is  the  frequency  of  the  downmixer.  The 
downmlxer  circuit  permits  decreasing  of  the  effective  freq¬ 
uency  shift  to  values  within  the  1  to  10  MHz  range.  The  true 
frequency  data  Is  then  averaged  over  the  samples  taken  to 
obtain  a  mean  value  and  a  sample  standard  deviation.  Using 
this  mean  value,  for  each  channel,  all  samples  that  are 
within  2  v,  two  standard  deviations,  of  the  mean  are  saved 
and  the  rest  are  discarded.  This  procedure  accounts  for  a- 
bout  5%  reduction  of  data  and  will  insure  samples  that  are 
too  low  or  high  will  not  bias  the  final  results.  A  new  mean 
and  standard  deviation  are  calculated  for  each  channel.  Next 
the  probe  volume  dimensions  and  parameters  are  calculated 


with  primary  emphasis  on  the  fringe  spacing  equations 

<  =  X  /  [2  #  sln<e/2>] 

where  6  Is  the  fringe  spacing  In  microns,  X  Is  the  beam 
wavelength  In  microns  and  ©  Is  the  angle  of  Intersection 
of  the  beams.  Knowing  the  fringe  spacing,  the  velocity  can 
be  determined  using  the  relations 

u  -  fd  #  * 

where  u  is  the  velocity  of  the  particle  in  meters  per  second 
and  is  frequency  of  the  particle  measured  by  the  count¬ 
ers  in  MHz. 

Once  the  velocity  components  of  the  particle  are 
determined  streamwlse  and  In  the  transverse  directions,  the 
computer  stores  the  results  in  a  data  table  for  further 
analysis.  The  probe  position  Is  also  saved  in  this  data 
file.  The  entire  process  is  now  repeated  at  different 
locations  of  Interest  In  the  flow. 

C.  LABORATORY  APPARATUS 

Besides  the  Laser  Doppler  veloclmeter  and  schlieren 
systems  described  earlier,  additional  instrumentation  was 
used  to  obtain  the  wind  tunnel  pressure  and  temperature 
data.  A  24  port  scan  1  valve  pressure  transducer  system  was 
set  up  to  measure  the  total  and  static  pressures  at  various 
points  Inside  the  operating  tunnel.  This  scanlvalve  system 
was  driven  by  a  digital  to  analog  Interface,  built  on  the 


premises,  which  allowed  computer  control  of  the  "stepping" 
and  "homing"  of  the  scanlvalve. 

Using  the  Digital  to  Analog  interface  on  the  MlcroVAX 
computer,  these  signals  to  step  the  scanlvalve  system  were 
transferred  from  a  computer  program  to  the  physical  unit  for 
action.  Once  the  stepper  motor  was  sent  the  appropriate  sig¬ 
nal  to  move,  the  Analog  to  Digital  Interface  board  on  the 
computer  sampled  the  pressure  transducer  signals  and  loaded 
the  data  Into  the  computer  for  conversion  Into  pressure  data 
and  later  analysis. 

The  atmospheric  pressure  was  recorded  from  a  mercury 
barometer  with  micrometer  scale  mounted  on  the  laboratory 
wall  near  the  wind  tunnel.  The  total  temperature  was  record¬ 
ed  on  a  mercury  thermometer  mounted  on  the  barometer.  Both 
of  these  measurements  will  soon  be  replaced  with  Paro- 
sclentlfic  pressure  sensors  and  thermocouples  that  will 
have  a  direct  link  with  the  computer. 

A  TSI  Six-Jet  Atomizer,  Model  9306,  was  used  to  provide 
particle  seeding  for  the  LDV  system. 

Figure  2.9  shows  the  general  laboratory  setup  and 
location  of  the  data  processing  and  ancillary  equipment. 

D.  DATA  PROCESSING 

1 .  Introduction 

The  Laser  Doppler  velocimeter  yields  a  large  amount 
of  raw  data  during  Its  operation.  This  data  must  be  stored 
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In  the  computer  In  a  logical  sequence  In  order  to  properly 


calculate  the  velocities  and  all  other  associated  data. 

2.  Computer  Programs 

The  following  discussion  pertaining  to  the  func¬ 
tioning  of  each  of  the  computer  programs  that  was  written  to 
accomplish  the  task  of  managing  the  vast  amount  of  data  In 
the  most  efficient  and  timely  manner. 
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a.  Master  Program 


MASTER  controls  the  overall  flow  of  data  during 
each  of  the  subroutine  calls,  acquires  all  the  data  and 
reduces  It  to  a  usable  form. 

First,  the  wind  tunnel  conditions  are  sampled  by 
the  use  of  various  subroutines:  DA11J1LAB,  SETUP,  GETDATA, 
DETACH  and  OUTPUTD.  The  traversing  mechanism  13  then  posi¬ 
tioned  with  subroutine  HWTRVT  to  move  the  probe  volume  of 
the  LDV  to  the  proper  location  in  the  flow.  Next,  the  LDV 
data  is  acquired  with  subroutine  HWWA.  The  samples  are 
taken,  velocities  calculated  using  subroutine  VELOCITY  and  a 
statistical  mean  and  standard  deviation  are  calculated.  The 
data  Is  then  reduced  with  subroutine  REDUCE_DATA  to  preclude 
the  use  of  any  velocities  greater  or  less  than  two  standard 
deviations  away  from  the  mean.  Each  of  the  two  velocity  com¬ 
ponents  are  then  recorded  along  with  the  position  In  the 
flow.  The  velocity  profiles  are  obtained  from  the  data  using 
post-acquisition  processing  programs.  An  option  Is  given 
during  the  acquisition  process  to  view  a  histogram  of  the 
resulting  velocities  to  check  on  the  distribution  of  data 
samples  with  subroutine  HISTOGRAM.  A  sample  histogram  from  a 
test  run  Is  shown  in  Figure  2.10. 
b.  Scan l valve  Programs 

First  a  signal  Is  sent  to  the  scanlvalve  to 
"home"  to  port  0  for  the  first  pressure  reading  This  Is 
accomplished  by  using  the  subroutine  DA11J1LAB.  Next,  a  call 
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Figure  2.10  Sample  Histogram 


to  subroutine  SETUP  Is  made  to  connect  the  computer's  In¬ 
ternal  timing  c  I  ock  and  Analog  to  Digital  interface  board  to 
program  control.  Also,  the  sample  acquisition  rate  and  num¬ 
ber  of  samples  are  interactively  entered. 

Now,  the  computer  is  ready  to  acquire  data  from 
the  scanlvalve  system.  Since  there  are  24  ports  to  sample  on 
the  scanlvalve,  a  loop  is  entered  in  which  samples  are 
taken,  scanlvalve  stepped  and  data  reduced  during  each  loop. 
Subroutines  GETDATA  and  DA11J1LAB  acquire  the  data  and  step 
the  scanivalve  respectively.  DA11J1LAB  has  an  adjustable 
time  delay  to  give  time  for  the  pressure  in  each  port  to 
equalize  after  scanivalve  stepping. 


After  all  24  porta  are  sampled,  the  mean  and 
standard  deviation  voltages  are  calculated  for  each  port. 
The  calibration  pressure  attached  to  port  0  and  the  atmos¬ 
pheric  pressure  at  port  1  are  now  entered  Interactively  In 
consistent  units;  inches  of  mercury  In  this  case.  With  these 
two  pressures  entered,  a  calibration  line  Is  calculated  and 
the  remaining  voltages  from  the  scanivalve  system  are  con¬ 
verted  to  pressures.  Using  port  2,  settling  chamber  stagna¬ 
tion  pressure,  and  port  7,  test  section  static  pressure,  a 
tunnel  Mach  number  is  calculated.  See  section  III  A  for  the 
equations  used  in  these  calculations.  Next,  the  total 
temperature  is  entered  and  the  free  stream  velocity  Is 
calculated. 

For  use  In  labelling  a  program  generated  plot¬ 
ting  file  for  post-acquisition  analysis,  the  angle  of  attack 
of  the  test  section  airfoil  and  current  time  are  entered. 

The  last  part  of  the  scanivalve  data  acquisition 
programs  close  the  loop  on  the  analog  to  digital  Interface 
by  detaching  the  clock  and  driver.  A  diagnostic  output  file 
Is  also  generated  for  post-acquisition  analysis, 
c.  LDV  Data  Acquisition  Program 

Both  the  transmitting  and  receiving  optics  of 
the  Laser  Doppler  Velocimeter  are  mounted  on  three  axis 
motorized  traverses  that  are  slaved  together  to  move  the 
probe  volume  to  any  point  of  Interest  in  the  flow.  Using  a 
pre-existing  software  driver  for  the  Taskmaster  traverse 
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control  unit,  the  optics  could  be  moved  at  will  from  the 
computer  workstation. 

Although  a  pre-existing  software  driver  was 
available  for  the  NASA-LDV  Computer  Interface,  significant 
changes  had  to  be  made  to  the  code  for  it  to  work  on  the 
present  system.  The  main  Laser  Doppler  Velocimetry 
acquisition  subroutine  "DAMUX"  was  responsible  for  loading 
the  incoming  data  from  the  multiplexer  into  data  arrays. 
Once  the  raw  frequency  data  is  in  arrays,  it  has  to  be 
decoded  to  correspond  to  usable  numbers. 

The  time  word  output,  as  14  bits,  by  the  counter 
has  a  mantissa  and  exponent  format:  with  a  10  bit  mantissa 
and  a  4  bit  exponent.  The  computer  reads  it  as  a  16  bit  word 
and  adds  two  bits  that  are  not  used  for  data  handling.  These 
two  bits  are  only  used  during  direct  computer  links.  This  16 
bit  word  is  output  as  a  complementary  word  that  is  reversed; 
that  is,  the  zeros  are  output  as  ones  and  the  ones  as  zeros 
then  the  entire  word  is  reversed  (l.e.  the  least  significant 
bit  which  is  on  the  right  most  side  bit  is  moved  to  the  left 
most  bit  in  the  word).  First,  the  word  is  made  into  its  com¬ 
plement  then  both  the  mantissa  and  the  exponent  are  reversed 
prior  to  decoding.  The  exponent  is  preset  from  the  counter 
front  panel  switch.  The  following  is  the  standard  16  bit 
word  composition  and  what  actually  results  from  the  output: 


Correct  output. 

0111010001  0111  01  =  465x27  =  59,520 

Mantissa  Exp  NU 


Actual  output 


01  0001  0111010001 


NU  Exp  Mantissa 


Translated  by  the  computer 
0100010111  0100  01  =  279x24  =  4,464 


First  step:  Take  its  complement 
10  1110  1000101110 
NU  Exp  Mantissa 


Next:  reverse  the  order 
0111010001  0111  01  =  59,520 

Mantissa  Exp  NU 

where  NU  stands  for  Not  Used.  In  order  to  fill  the  sixteen 
bit  word,  the  14th  and  15th  bit  are  added  to  control  data 
flow  but  have  no  significance  to  the  value  of  the  data.  From 
this  last  16  bit  word  is  calculated  the  Doppler  shift  freq¬ 
uency.  The  output  from  the  NASA-LDV  Computer  Interface  is 
actually  representative  of  the  time  elapsed  for  the  particle 
to  traverse  8  fringes  of  the  probe  volume.  The  Doppler 
period  is  expressed  by  the  following  equation: 


td  =  1/32  (D9D8...D0)  x  2®^  ns 

where  D9..D0  represents  the  binary  mantissa  and  EXP  is 
the  preselected  exponent  set  on  the  counter.  From  the 
period,  the  Doppler  frequency  is  calculated  using  the 
following  equation: 

fd  =  1000/td  MHz 

This  Doppler  frequency  is  then  be  used  to  calculate  the 
velocity  of  the  traversing  particle.  Since  frequency 
shifting  is  used  to  discriminate  the  actual  direction  of  the 
particle  through  the  probe  volume,  the  Doppler  frequency 
must  be  adjusted  to  account  for  the  Bragg  and  mixing 
frequencies.  See  equation  on  page  34  for  this  relationship. 
Once  the  actual  frequency  of  the  particle  is  determined,  the 
fringe  spacing  is  calculated  by  the  first  equation  on  page 
35  and  the  velocity  is  found  using  the  second  equation  on 
the  same  page. 
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d.  LDV  Support  Programs 

Three  subroutines  are  used  to  provide  the 
supporting  functions  to  the  LDV  data  acquisition.  These  are 
VELOCITY,  HISTOGRAM  and  REDUCE_DATA . 

VELOCITY  is  responsible  for  calculating  the 
fringe  spacing,  extracting  the  particle  frequency  from  the 
signal  and  then  calculating  the  two  velocities. 

HISTOGRAM  takes  the  velocity  data  and  plots 
a  histogram  of  the  percentage  of  velocities  that  fall  within 
a  given  range.  The  x  axis  of  the  histogram  Is  divided 
equally  Into  8  segments  from  the  minimum  to  maximum  velocity 
values.  The  option  to  plot  a  velocity  histogram  Is  given 
after  each  LDV  acquisition  to  spot  check  the  distribution  of 
the  data.  A  nearly  Gaussian  or  normal  distribution  Is 
expected  if  the  system  Is  operating  properly.  This  will  also 
show  the  effects  of  noise  as  evidenced  by  the  broad  band  of 
the  histogram  values. 

REDUCE_DATA  reduces  the  velocity  data  to 
only  those  values  falling  within  two  standard  deviations  of 
the  mean.  This  results  In  a  4.5%  rejection  of  data  which 
will  reduce  the  effect  of  high  frequency  noise  on  the 
counters  and  tend  to  reduce  the  effect  of  very  large 
particles  moving  slowly  through  the  probe  volume. 


A.  TUNNEL  CONDITIONS 


In  order  to  obtain  good  correlation  between  the  schlie- 
ren  photographs  and  the  Laser  Doppler  Veloclmeter  obtained 
data,  the  operating  tunnel  conditions  would  have  to  be 
nearly  identical.  Two  sets  of  conditions  were  chosen:  one 
corresponding  to  an  approximate  Mach  number  of  0.3  and  the 
other  0.5.  Once  the  tunnel  was  running  at  operating 
conditions,  the  computer  program  MASTER  calculates  the  test 
section  Mach  number  and  free  stream  velocity  according  to 
the  following  relationships.  Assuming  lsentroplc  flow  in  the 
test  section,  the  section  Mach  number,  M,  is  calculated  from 
the  stagnation  pressure,  PQ,  and  the  static  pressure,  P, 
using  the  following  relationship: 

M  -  (aP0/P]<Tf_1/v>-l)x(2/[r-13>)1/2 

where  y  =  1.4.  The  resulting  calculated  Mach  numbers  for 

the  two  conditions  are  as  follows:  0.295  and  0.450.  From 
these  Mach  numbers,  the  corresponding  free  stream  velocities 
are  calculated. 

V  »  M  #  a 

where  a  is  the  sonic  velocity. 

a  -  (yRT>1/2 

where  R  is  the  universal  gas  constant  and  T  is  the  static 


temperature  In  the  test  section.  The  static  temperature  can 
be  calculated  from  the  total  temperature,  TQ,  measured  on 
the  wall  near  the  tunnel  by  the  next  relationship: 

T  =  Tq/U  +  [  (  7-1  >/2]#M2) 

B.  SCHLIEREN  PHOTOGRAPHS 
1 .  Introduction 

Once  the  wind  tunnel  was  running  at  the  first  of 
two  operating  conditions,  manometer  board  readings  were 
recorded  and  Mach  number  calculated  using  the  first 
equation  on  the  prelous  page.  The  schlleren  photographs  were 
then  taken  at  four  different  angles  of  attack  and  three 
different  strobe  Intensities  to  obtain  the  optimal 
photographs. 
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2.  Test  Conditions  and  Procedure 

Using  the  procedures  detailed  In  section  II  A,  the 
schlleren  photographs  were  taken  at  the  settings  and 
conditions  listed  In  Table  I,  below. 

Upon  completion  of  the  wind  tunnel  runs,  the  35mm 
film  was  developed  and  printed.  From  these  schlleren 
photographs,  the  conditions  resulting  In  the  best  depiction 
of  the  wake  profiles  were  chosen  to  represent  each  Mach 
number  and  angle  of  attack.  These  are  annotated  with  an 
asterisk,  #,  In  Table  I,  below. 


Using  these  photographs,  a  short  HP-41  CX  calculator 


program  was  written  to  relate  photographic  dimensions  to 
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TABLE  I. 

SCHLIEREN  PHOTOGRAPH  PARAMETERS 

Strobe 

£1 

a  (deg) 

Intensl tv  (Hz 

0.30 

0.0 

20 

0.30 

0.0 

70 

0.30 

0.0 

150 

0.30 

7.7 

20 

0.30 

7.7 

70 

0.30 

7.7 

150 

0.30 

15.5 

20 

0.30 

15.5 

70 

0.30 

15.5 

150 

0.30 

23.6 

20 

0.30 

23.6 

70 

0.30 

23.6 

150 

0.45 

0.0 

20 

0.45 

0.0 

70 

0.45 

0.0 

150 

0.45 

7.7 

20 

0.45 

7.7 

70 

0.45 

7.7 

150 

0.45 

15.5 

20 

0.45 

15.5 

70 

0.45 

15.5 

150 

i 

r< 

1 

0 

r, 
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real  teat  section  dimensions.  Once  the  real  dimensions  of 
the  wake  were  known,  the  Taskmaster  traverse  mechanism  was 
calibrated.  This  yielded  a  1:1  relationship  in  the  chordwlse 
direction  and  a  1:4  relationship  In  the  vertical  direction. 
This  means  that  the  Taskmaster  digital  control  panel  reads 
true  values  for  the  X  direction,  chordwlse,  and  four  times 
the  real  value  In  the  2  direction,  vertical.  Arbitrarily 
choosing  the  trailing  edge  of  the  airfoil  at  0°  a  as  an 
Initialization  point  for  the  Taskmaster,  a  table  of  LDV 
sampling  points  was  made  for  the  first  three  angles  of 
attack.  See  Table  II,  below,  for  a  sample  of  the  Taskmaster 
table  at  a  ■  7.7° 
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Points 

SMPitd  X  2 


Alpfta  *  ?.?  deg 


20  2.00 

7 

6.5 

6 

5.5 

5 

4.5 

4 

3.5 

3 

2.5 

2 

1.5 

1 

.5 

0 

-.5  -1 

-1.5  -2  -2.5 

20  1.50 

7 

6.5 

6 

5.5 

5 

4.5 

4 

3.5 

3 

2.5 

2 

1.5 

1 

.5 

0 

-.5  -1 

-1.5  -2  -2.5 

20  1.00 

7 

6.5 

6 

5.5 

5 

4.5 

4 

3.5 

3 

2.5 

2 

1.5 

1 

.5 

0 

-.5  -1 

-1.5  -2  -2.5 

20  .50 

7 

6.5 

6 

5.5 

5 

4.5 

4 

3.5 

3 

2.5 

2 

1.5 

1 

.5 

0 

-.5  -1 

-1.5  -2  -2.5 

20  .00 

7 

6.5 

6 

5.5 

5 

4.5 

4 

3.5 

3 

2.5 

2 

1.5 

1 

.5 

0 

-.5  -1 

-1.5  -2  -2.5 

19  -.50 

7 

6.5 

6 

5.5 

5 

4.5 

4 

3.5 

3 

2.5 

2 

1.5 

1 

.5 

0 

-.5 

-1.5  -2  -2.5 

17  -1.00 

7 

6.5 

6 

5.5 

5 

4.5 

4 

3.5 

3 

2.5 

2 

1.5 

1 

.5 

-1.5  -2  -2.5 

17  -1.50 

7 

6.5 

6 

5.5 

5 

4.5 

4 

3.5 

3 

2.5 

2 

1.5 

1 

.5 

-1.5  -2  -2.5 

11  -2.25 

7 

6.5 

6 

5.5 

5 

4.5 

4 

3.5 

3 

2.5 

-2.5 

12  -3.00 

7 

6.5 

6 

5.5 

5 

4.5 

4 

3.5 

3 

2.5 

2 

-2.5 

176  (Total  locations  sapleti) 


C.  LASER  DOPPLER  VELOCIMETRY  DATA  ACQUISITION 


The  relationship  between  the  Doppler  frequency  and 
the  velocity  of  the'partlcle  traversing  the  probe  volume  is 
given  by  the  second  equation  on  page  35.  The  particle 
velocity  is  directly  proportional  to  the  Doppler  fre<*iency 
by  a  factor  of  the  fringe  spacing.  By  the  first  equation  on 
page  35,  fringe  spacing  is  only  a  function  of  the 
transmitting  lens  angle  and  wavelength  of  the  light.  Since 
the  wavelength  of  the  light  is  accurately  known  to  within 
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0.05  microns,  or  0.01\,  an  accurate  lens  angle  Is  all  that 
Is  needed  for  calibration  [Ref.  41. 


The  following  relationship  was  used  to  calculate 
the  lens  angle: 

e  =  2  *  tan-1  (E#d/(2#  f>) 

where  E  is  the  beam  expansion  ratio,  1.0  in  this  case:  d  is 

the  incident  beam  spacing,  preset  at  13  mm  with  TSI  Model 

9113-13?  and  f  is  the  wavelength  of  the  light. 

The  lens  angle  in  this  case  was  not  accurate 

because  the  focal  length  of  the  focusing  lens  was  not 
known.  Therefore,  calibration  was  necessary  to  insure 
accurate  velocity  measurements  were  obtained. 

The  calibration  method  chosen  was  to  operate  the 
wind  tunnel  at  known  conditions,  calculate  the  free  stream 
velocity  from  the  scanlvalve  pressure  data  and  then  move 
the  probe  volume  via  the  Taskmaster  to  coincide  with  a 

position  In  which  the  free  stream  velocity  could  be 
measured.  From  this  procedure,  a  constant  could  be  applied 
to  the  velocity  equation  as  follows: 

u  =  fd  #  i  #  (Calibration  Constant) 

This  calibration  constant  will  compensate  for  the  estimated 
focusing  lens  focal  length. 

2.  Test  Conditions  and  Procedure 

Once  calibration  Is  completed,  the  tunnel  Is 
brought  to  the  test  conditions,  the  Taskmaster  is  moved  to 


the  first  sample  point  and  the  MASTER  computer  program  Is 
run  with  the  adjusted  fringe  spacing  to  reflect  the 
cal Ibrat Ion . 

At  each  sampling  point,  1000  samples  will  be  taken 
to  Insure  statistical  accuracy  of  the  mean  flow  calcula¬ 
tions.  The  scheduled  sampling  points  can  be  seen  In  Table 
II,  page  48.  After  each  Taskmaster  traverse,  the  oscillo¬ 
scope  Is  checked  to  insure  no  drastic  loss  of  the  signal 
occurred  during  the  movement.  If  the  Doppler  signal  de¬ 
creases,  slight  adjustments  to  the  receiving  optics  traverse 
may  be  needed  to  regain  the  signal  strength.  To  optimize  the 
data  received,  slight  adjustments  may  have  to  be  made,  also, 
to  the  selected  mixing  frequency,  mixing  frequency  ampli¬ 
tude,  counter  threshold  setting  or  high  and  low  pass  filter 
sett Ings. 

After  Table  II  data  points  are  sampled,  the  airfoil 
angle  of  attack  or  test  Mach  number  Is  changed  and  the  above 
procedure  Is  repeated. 


A.  SCHLIEREN 


The  first  thesis  goal  was  attained  early  enough  to 


provide  qualitative  Information  concerning  the  airflow 


around  and  in  the  wake  of  the  airfoil.  Schlleren  35mm 


photographs  were  produced  for  both  selected  mach  numbers. 


0.3  and  0.45,  and  for  four  different  angles  of  attack;  0°, 


7.7°,  15.5°  and  23.6°.  The  best  photographs  for  each 


Mach  number  and  the  first  three  angles  of  attack  are  dis¬ 


played  in  Appendix  A,  pages  55-57.  From  these  photographs. 


LDV  sampling  points  were  chosen  as  seen  from  Table  II  on 


page  48. 


Three  turbulent  areas  are  evident  on  all  the  schlleren 


photographs.  As  seen  in  Figure  A. 4  on  page  56,  the  airflow 


is  slightly  disturbed  upstream  of  the  airfoil  due  to  the 


presence  of  a  pitot-static  tube  extending  into  the  airflow 


Just  forward  of  the  airfoil.  This  tube  has  been  subsequently 


removed  to  clean  up  the  flow.  At  the  top  of  the  same  photo¬ 


graph,  a  wedge  shaped  disturbance  can  be  seen.  This  and  the 


ower  floor  turbulent  layer  are  results  of  the  transonic 


wall  slottings  found  on  both  the  roof  and  the  floor  of  the 


test  section.  Although  a  concern  for  later  unsteady  flow 


studies,  these  two  areas  do  not  Interfere  with  the  wake 


profiles  resulting  from  a  steady  flow  analysis. 
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B.  LASER  DOPPLER  VELOCIMETER 

Due  to  equipment  problems  and  time  constraints,  the  LDV 
data  acquisition  part  of  this  thesis  was  not  able  to  be 
accomplished.  All  systems  are  now  operational,  to  include 
computer  programs,  however,  there  is  inadequate  time  to 
acquire  data  and  analy2e  the  velocity  profiles  obtained. 

Calibration  runs  were  made  with  a  limited  degree  of 
success  when  the  wind  tunnel  was  run  under  Mach  0.3  and  the 
Laser  was  not  powered  over  2  watts.  The  problem  with  the 
Laser  is  easily  repaired,  however,  time  consuming.  The 
upper  limit  on  the  wind  tunnel  Mach  number  is  caused  by  the 
evaporation  of  the  seeding  particles  in  high  speed  flow. 
This,  too,  can  be  easily  overcome  by  the  proper  selection 
of  seeding  particles  and  related  atomizer. 
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A.  CONCLUSIONS 


The  thesis  provided  Invaluable  assistance  to  the  dynamic 
stall  project  at  the  Fluid  Mechanics  Laboratory.  The  schlie- 
ren  photography  was  able  to  provide  the  qualitative  data  on 
the  steady  state  flow  field.  With  this  accomplished,  the 
schlieren  system  in  this  laboratory  was  able  to  be  easily 
tailored  to  the  next  step:  unsteady  flow  analysis. 

The  LDV  work  provided  the  laboratory  with  the  essential 
computer  programs,  hardware  Interfaces  and  drivers,  optics 
systems  and  auxiliary  required  equipment  to  enable  the  pro¬ 
ject  to  proceed  with  very  little  additional  work  to  achieve 
the  steady  state  velocity  profiles.  The  ground  work  has  been 
accomplished  to  enable  the  project  to  continue  forward  with 
the  systems  on  hand  to  study  the  effects  of  dynamic  stall  in 
both  the  steady  and  unsteady  environments. 

B.  RECOMMENDATIONS 

The  schlieren  system  is  fully  functional,  however,  the 
Laser  Doppler  Veloclmeter  system  still  requires  a  little 
additional  effort  to  make  it  productive.  The  Laser  doughnut 
condition  needs  to  be  addressed  to  enable  higher  powers  to 


be  utilized.  The  types  of  seeding  available  to  the  system 
needs  to  be  addressed  to  insure  operability  of  the  tunnel 
and  data  acquisition  systems  above  Mach  0.3.  With  these 
minor  fixes  and  a  little  fine  adjustments  of  the  signal 
processing  equipment,  the  LDV  Is  ready  for  action. 
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Figure  A. 3 


Mach  Number 


Figure  A. 4  Mach  Number 
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APPENDIX  Bs  COMPUTER  PROGRAMS 


PROGRAM  MASTER 

- - - - —  — - - 

C  Most e r  Progrom  written  by  Moj  Richord  R  Rylss,  NPS  Student  28  Apr  87 
C  Lost  Update  6  Ssp  87 

C 

C  This  program  will  gather  data  from  the  scanivolve  system,  move  the 
C  Taskmaster  and  acquire  LDV  velocity  dota. 

C  This  program  should  be  used  with  subroutine  SETUP.  GETDATA,  DETACH.  OUTPUTD 
C  HWTRVT.  HWWA,  HISTOGRAM  and  VELOCITY. 

C - - - 

INCLUDE  'SYS$L1BRARY:L10SET . FOR’ 

C - — - - - — - - - - 

C  Declare  variables 

C - - - - 

INTEGER  D£V_ID.CLK_ID 

INTEGER  NCHANL.HISTO.BLUE_GREEN 

INTEGER  NPT,NPTS,NSAU»LES. COUNT 

INTEGER  I. J.NSTEPS.IER. ANSWER. EH>. TALLY. TIME 

REAL  VOLTAGES( 5808) . U_U( 208 ) , V_V ( 208 ) , U_V( 288 ) . OELX . OELZ 

REAL»8  sum.amean,x.xx.y.yy.sto_oev.sto_dev  V.STD_DEV  U.SUM_V.SUM  U 

REAL'S  SUM  VERT.SIALHOR I Z 

REAL  RAT£.SCANI_VOLT(8:23.2),SCANI_PRESS(8:23).VERT  VEL(8208) 

REAL  CALI0_PRESS. SLOPE. INTERCEPT. HORIZ  VEL ( 8208 ).F I NAL_V. FINAL  U 
REAL  X_POSITION.*J>OSITION.M£AN(2).M.FREE_STREAM. TO. ATMOS  PRESS 

REAL  X.LOC ( 288)  .  Z_LOC( 288)  ,U(288) ,  V(288)  ,  ALPHA. ZERO. CORR_ FACTOR 

REAL'S  UU.W.UV 
CHARACTER.  1  ANS.VOLTjtMS 
LOGICAL  HOME 

C - - - - - - - 

C  Sconivalve  data  acquisition 

C - — - — — - - - 

0PEN(UNIT-99.  FILE-  RESULTS.  OAT ’.STATUS- ’NEW’) 

TALLY  -  8 
HOME  -  .TRUE. 

NS  T EPS  -  8 

CALL  0A11J1LAB  (HOME.MSTEPS. I  EH) 

IF  (IER.NE.1)  GOTO  178 
ANSWER  -  8 

WRITE  (•.')  ’Skip  to  Taskmaster?  (y/n)  [DEF  -  No]’ 

READ  (*.1)  ANS 

IF  (ANS.EQ. ’Y’ -OR.AMS-EO. ’y’)  GOTO  98 

WRITE  ('.•)  ’Print  the  sconivalve  diagnostic  data  to  a  file?  (Y/N)’, 

*  ’  [DEF  -  No]’ 

READ  («.1)  VOLT.ANS 

CALL  SETUP  (CLK_IO.DCV_ID.RATE.NCHAML.NPT.NPTS) 

18  CONTINUE 

HOME  -  .FALSE. 

NS  T  EPS  -  1 
DO  58  J  -  8.  23 
WRITE  (• .2)  J 

IF  (J.NE.0)  CALL  0A1 1 J 1  LAS  (HOME.MSTEPS . I ER) 

IF  (IER.NE.1)  GOTO  158 
28  CONTINUE 

CALL  GETDATA  (DEV_ID.NPT.»»TS. VOLTAGES) 

C  Calculate  eean  and  standard  deviation 

. . . . . 

SUM  -  8.0 

DO  38  I  -  2.  NPT 

SUM  -  SUM  +  VOLTAGES ( I ) 


r 


58 


30  C0NT1NUC 

AMCAN  -  SUM  /  FLOAT (NPT-1) 

SUI  -  0.D0 
DO  40  I  -  2.  NPT 

X  -  VOLTAGES ( I )  -  AMCAN 
Y  -  X  •  X  •  0  .0/FLOAT(W>T-1)) 

46  SUM  -  Y  4  SUM 

STO.OEV  -  SORT (SUM) 

IF  (VOLT  ANS.EQ. 'y  ’ . OR . VOLT_ANS . EQ. 'Y'  )  CALL  OUTPUTO  (RATE. NPT, 
t  NCHANL.AMEAN.STD.DEV. J, VOLTAGES) 

SCAN  I _VOlT ( J .  1 )  -  AMCAN 
50  SCANI_V0LT(J.2)  -  STO.OEV 

CALL  OETACM  (CLK_IO,OEV_ID) 

IF(VOLT_ANS.EQ.  *Y* .  OR .  VOLT_ANS .  EO .  'y')  THEN 

OP EN(UNIT-1 5. FILE- 'SC ANI_VOLT.DAT1 . STATUS-' NEW ' ) 

00  60  1  -  0.  25 

60  HWtlTE  (15.5)  I.  SCANI_V0LT(I.1).  SCANI_VOLT( 1 .2) 

END  IF 

C********************** 

C  Calculate  ecanivalve  calibration  lino 
C  and  convart  all  voltogee  to  preeeure 

. . . . . 

WRITE  (*,*)  'Enter  Scani val ve  Col ibrot I  on  Preeeure:’ 

READ  !*.*)  CALIB_PRESS 

WRITE (•.*)'Enter  Calibration  Preeeure  Correction  Foctor:' 

READ( • , • )CORR_F ACTOR 

CALIB.PRESS  -  CALIB.PRESS  +  CORR_ FACTOR 
WRITE(* ,*)' Enter  atiaoapheric  preeeure:' 

READ( • , • )  ATMOS_PRESS 

SLOPE  -  (ATMOS_PRESS-CALIB_PRESS)/(SCANI_VOLT  (1.1  )-SCANI_VOLT(0, 1 ) ) 

INTERCEPT  -  CALIB_PRESS  -  SCANI_VOLT(0. 1 )»SLOPE 

WRITE(99. •)  CALIB.PRESS' .CALIB.PRESS 

WRITE! 99. •)' CORR_F ACTOR ' ,CORR_F  ACTOR 

WRITE! 99, •)  ’ ATMOS_PRESS ' .ATMOS.PRESS 

OP  EN(UNIT-10,nLE-' SCAN  IPORT.DAT' .  STATUS- '  NEW  *  ) 

DO  70  I  -  0.  25 

SCAN I .PRESS  (I)  -  SCAN I _ VOLT (1.1)  •  SLOPE  ♦  INTERCEPT 
70  WRITE( 16 . • ) I ,  SCANI_PRESS(I) 

CLOSE( 16) 

IF  (VOLT.ANS.EQ.'Y' .OR.VOLT_AM5.EQ.'y’)  THEN 
DO  80  I  -  0.  23 

00  WRITE  (15.4)  I ,SCANI_PRESS( I ) ,SCANI_VOLT ( 1 . 2) •SLOPE+INTERCEPT 

EfOIF 
CLOSE! 15) 

WRITE! • .  •  ) '  ' 

WRITE!* , •) 'Equot ion  of  the  calibration  line  •  * 

WR1TE(*.*)'Y  -'.elope,'  X  . intercept 
WRITE(99. •) *Y  -'.elope.'  X  intercept 

write!*. •)'  ' 


C  Calculate  tunnel  Mach  and  velocity 

C  •  •  •  •  •  . . . 

IF(SCANI_PRESS(7).EQ  0.0)  SCANI_PRESS(7)  -  1.0 
IF(SCANI_PRESS!2) . LT  SCAN 1 .PRESS ( 7 ) )  GOTO  85 
M  -  SQRT(((SCANI_PRESS(2)/SCANI_PRESS(7))**(2/7.)-l  )*5.) 
85  WRITE! *.*) 'Tunnel  Mach  Number  -'.M 
WRITE(99,*) 'Tunnel  Mach  Nieeber  -' ,M 
VfRITCt  •  •  )  • 

WRITEl* .*) 'Enter  Total  Temperature  (C) ’ 

READ(* . •)  TO 

WRITE(99.*) 'TOTAL  TE>0>'  ,  TO 
WRITE(* , •) '  ' 

TO  -  TO  ♦  273.15 

FREE_STREAM  -  M  •  SQRT(401 . 52*T0/(  1 .  4(M**2/5 .  ) ) ) 

WRITE!* .•) 'Free  Stream  Velocity  (mpe)  -' . FREE.STREAM 
WRITE! 99. «)'Free  Stream  Velocity  (mpe)  -'.FREE.STREAM 
IF(FREE_STREAM.EQ.0.0)  FREE.STREAM  -  IE-10 
WRITE!*. •) '  ’ 


00000-000  OOO  —  40000 


WRITE(» . •) ' Enter  angle  of  attack  (dog):' 

RCAO(* . •)  ALPHA 
WRITE(99.  •)  'ALPHA'  .  ALPHA 

WRITE(«.«)’Enter  Curront  Tine:  (4  digit  intogor)' 

Reod(* . •)  TIME 
WRITE(99.«) 'TIME*. TIME 

WRITE  («.•)  'Skip  to  ond?  (y/n)  [DEF  -  No]' 

READ  (*.1)  ANS 

IF  (ANS .  EO .  ' y  ’  . OR .  ANS .  EO .  ‘ Y*  )  GOTO  1M 

Taokmooter  aover  ond  LDV  dato  acquisition  | 

WRITE  (*.*)  ‘Which  color  it  connoctod  to  tho  loft  (|3)  channol?'. 
t  '  (81 ue-1 , Green-9)' 

READ  (*.9)  BLUE.GREEN 

WRITE(99. •) 'BLUE_GREEN' ,8LUE_ GREEN 

I F ( FR£E_STR£AM . EO . 6 . 6)  FREE_STREAM  -  1 . E-« 

OPEN  (UNIT-25. FILE- ’FINAL. OAT’ .STATUS- 'NEW') 

WRITE  (25, •  )  *  Taekmoeter  (Final. Dot)’ 

WtlTE  (25. •)  ’  poo  (in)  velocity/etd  dov  velocity/etd  dov’ 

WRITE  (25,16) 

WRITE  (•.•)  'Skip  to  LDV?  (y/n)  [Def  -  NO]’ 

READ  (*.1)  ANS 

IF  (ANS . EQ . ’ Y ’ . OR . ANS . EQ . ’ y ‘ )  GOTO  119 
IF(TALLY.EO.e)  THEN 

WRITE(* . •) ’ Entor  X-Po«ition  for  calibration' 

READ(«.»)  X.POSITION 

WRITE(» .»)’ Entor  Z-Position  for  calibration* 

READ(«.*)  Z_POSITION 
ENOIF 

CALL  HWTRVT  (DELX.DELZ) 

WRITE  (*.«)  'Movo  tho  Taokioaotor  again  prior  to  LDV  Acquioition?  (Y/N)' 
READ  («,1)  ANS 


Rocord  pooition  of  LDV  at  all  tiwao 


IF (ANS . EQ. ’ Y' .OR. ANS . EQ . ’ y ’ )  THEN 
X.POSITION  -  X.POSITION  ♦  OELX 

Z.POSITION  -  Z_POSITION  ♦  DELZ 
GOTO  199 

EfCIF 


Acquiro  LDV  data 


19  CALL  HWWA  (NSAAFLES. ANSWER ,BLUE_CREEN) 

Increment  pooition  and  rocord  it 
Road  velocity  dato,  stoo  it  and  find 
otandord  deviation 


X.POSITION  -  X_POSITION  +  DELX 
Z.POSITION  -  Z_POSITION  +  DELZ 
SUBVERT  -  9.09 
SUH.HORIZ  -  9.09 

OPEN  (UNIT-35. FILE-' VELOCITY. DAT' , STATUS- ’ OLD ’ ) 

DO  129  I  -  1.  NSAJ*>LES 

IF  (8LUE.GREEN . EO . 9)  READ(35.7)  VERT_VEL( I ) ,HORIZ_VEL( I ) 
IF  (BLUE.GREEN. EQ. 1 )  READ(J5.7)  HORIZ_VEL( 1 ) , VERT_VEL( I ) 
SUM.VERT  -  SUM.VERT  ♦  VERT_VEL(I) 

129  SUM.HORIZ  -  SUM.HORIZ  ♦  HORIZ_VEL(I) 

CLOSE (35) 

FINAL.V  -  SUBVERT  /  FLOAT (NSAJyPLES) 

FINALJJ  -  SUM_HORIZ  /  FLOAT (NSAfcFLES) 

SUMJf  -  9.09 
SUM_U  -  9.D9 
UU  -  9.09 


ooo  o  o  o 


00  IN  I  •  t,  NSAIPl.CS 

x  -  vorr_vtL(i)  -  final_v 

XX  -  HORI Z_VEL(  1 )  -  FINAL.U 
Y  -  X  •  X  *  ( 1 .  #/FL0AT (NSAAFLES— 1 ) ) 

YY  -  XX  •  XX  •  (1  .6/FL0AT(NSAAPLES-1)) 

UU  -  UU  ♦  ((HORIZ_VEL( I )**2)-(FINAL_U**2) ) 

W  «  W  ♦  ??VERT_VEL(I)**2)-(FINAL_V**2)) 

UV  -  UV  ♦  ? ?HORIZ_VEL( I )*VERT_VEL( I ) )-(FINAL_U»FINAI _ V) ) 

SUM.V  -  Y  ♦  SUM.V 
136  SUM.U  -  YY  ♦  SUM.U 

UU  -  UU  /  FLOAT (NSAAPLES) 

W  «  W  /  FLOAT  (NSAJPLES) 

UV  -  UV  /  FLOAT (NSAAPLES) 

STO_DEV_V  -  SORT  (SUM.V) 

STO_OEV_U  -  SORT  ?SUM_U) 

WRITE?*.*)  '  1 

WRITE  (• . •) 'Ve loe i t ies  and  std  d«v» 

WRITE  (•.*)'U  di r'  ,FINAL_U.STO_DEV_U 
WRITE  (• .  •)  ’V  dir' ,FINAL_V,STD_DEV_V 
WRITE  (99, •) ‘Vel oe i t i es  and  «td  devs  » ’ 

WRITE  (99 ,  •)  ’U  d i r ’ . F I NAL.U . STD.DEV.U 
WRITE  (99,«)'V  d i r ' , F I NAL_V . STD.DEV.V 
WRITE  (*.*)’  ‘ 

WRITE?* . •) ' F i rat  cut  RMS  valuaa' 

WRITE?* . •)  '  U-RMS .  V-RMS  ‘  .UU.W 
WRITE?*. *)'UV-RMS’.UV 
WRITE?99.*)'Fi rat  cut  RMS  valuaa' 

WRITE?99,*J  ’U-RMS.V-RMS'  .UU.VV 
WRITE(99. •) ’UV-RMS' .UV 


El  i  Mi  not#  data  outaida  tare  aigwa  from  Mean 


IF  (BLUE.GREEN.EQ.1)  THEN 
MEAN?1)  -  FINAL_U 
MEAN? 2)  -  FINAU.V 

CALL  REDUCCJJATA  (NSAAFLES , COUNT , ENO , MEAN . STD_DEV_U , STD.DEV.V . 
k  BLUE.GREEN) 

F I NAL_U  •  MEAN(1) 

FINAL_V  *>  MEAN? 2) 

WRITE  (99.*) ‘Reduced  velocities  and  std  devs  ■’ 

WRITE  (99.*) 'U  dir' , F I NAL_U , STD_DEV_U 
WRITE  ?99. •) ’V  dir' . FINAL_V , STD_DEV_V 

ELSE 

MEAN?2)  -  FINALS 
MEAN? 1 )  -  FINAL. V 

CALL  REDUCE.OATA  (NSAfcPLES. COUNT, END, MEAN. STD_DEV_V . STO_OEV_U . 
ft  BLUE.CREEN) 

FINAL.U  -  MEAN? 2) 

FINAL. V  -  MEAN? 1 ) 

WRITE  (99, •) 'Reduced  velocities  and  etd  devs  »’ 

WRITE  ?99,*)'U  dir' . FINAL.U, STO.OEV.U 
WRITE  (99. •) 'V  di r' .FINAL_V.STD_OEV_V 

ENO  IF 


Calculate  RMS  values 


OPEN (UNI T»43, FI LE“’ VELOCITY . OAT  * . STATUS* ' OLD  * ) 

UU  -  6.06 
W  -  6.06 
UV  -  6.06 
00  216  I  -  1.  ENO 

IF  (BLUE.CREEN.CQ.6)  REAO?43,7)  VERT_VEL( I ) .HORIZ.VEL? I ) 
IF  (BLUE.CREEN. EO. 1)  REAO?«3.7)  HORIZ_VEL( I ) .VERT.VEL? I ) 
UU  -  UU  4  ( ?HORIZ_VEL( I ) **2)-( FINAL.U* *2) ) 

W  -  W  ♦  ??VERT_VEL(I)**2)-(F1NAL_V.*2)) 

216  UV  -  UV  *  (?M0RIZ_VEL(I)*VERT_VEL(1))-(FINAL_U*FINAL_V)) 
CL0SE(43) 

UU  -  UU  /  FLOAT (EM)) 
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UV  -  UV  /  FLOAT  (DC) 

WRITEf* . • V  ‘ 

WRITE! *.*) 'Rsducad  RMS  voluss' 

WRITE(«.«)'U-RMS,V-RMS‘ ,UU,  W 
WRITE! •.•j'UV-RMS*  ,UV 
WRITE(99. •)  ' U-RMS ,  V-RMS 1  ,UU,  W 
MITE(99 ,  •) ‘UV-RMS’ ,UV 

WRITE  (25.8)  X_POSITION.Z_POSITION.FIMAL_U.STO_OEV_U.FIHAL_V.STD  DEV  V. 
k  UU.W.UV 

TALLY  -  TALLY  ♦  1 


Options  to  visw  histogram  and  taka  doto  again 
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WRITE(«.«)'  ‘ 

MITE  (*.«)  ‘Do  you  want  a  Histogram?  (1/0)' 

READ  (•.•)  HIST0 

IF  (HISTO.EO. 1)  CALL  HISTOGRAM(COUNT.E».BLUE_GREEN) 
WRITE  ( • , 6) 

READ  (*.1)  ANS 

IF  (ANS. EO. 'n' .OR.ANS.EO. 'N' )  GOTO  140 
GOTO  100 
CLOSE  (25) 

CLOSE(99) 


Convsrt  output  data  into  propsr  format  for 
"PLOTCV"  Vsctor  Plot  of  vslocity  profilss 
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0PEN( UN I T-45 . F I LE- ‘ F I HAL . OAT ' . ST ATUS- ‘ OLD ‘ ) 

READ (45. 14) 

OP EN( UN IT-46. FILE- ‘CRA. DAT1 .STATUS- ‘NEW') 

WRITE(40.13)  M. FREE.STREAM, ALPHA, TIME. 11.10,10 
ZERO  -  00 
DO  150  1-1. TALLY 

REA0(45. 1 1 )X_L0C( I ) ,Z_L0C( I),U(I),V(I) .U_U( I ) . V_V( I ) ,U_V( I ) 

X_L0C ( I )-X_L0C( 1)43.0 
Z_L0C( I )-(Z_L0C( I )+2 . 5) • . 24 
U(I)  «  U(I)  /  FREE.STREAM 
V(I)  -  V(I)  /  FREE.STREAM 
U_U(I)  -  U_U(I)  /  FREE_STREAM 
V.VI  -  V_v(l)  /  FREE.STREAM 
U_V(I)  -  U_V(I)  /  FREE.STREAM 
WRITE(40,12)I, X_L0C( I ) . Z_L0C( I ).U(I),V(I) ,ZERO.U_U( I ) . V_V( I ) . ZERO. 
U_V( I ) . 7 


CLOSE (45) 
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CLOSE(40^ 
WRITE  (• .5) 


.ZERO 


Format  statsmsnts 


FORMAT 

FORMAT 

FORMAT 


FORMAT 


FORMAT 


(A1 ) 

(ix. 'Seanivalvs  f  *,I2) 

(IX. 'Port  #* . I2.3X, 'Avsrags  Voltogs  -‘.F12.4.3X, 

‘Standard  Dsviotion  -’.F12.4) 

(IX. ‘Port  #‘ . 12.3X, ‘Prsssurs  -\F12.4.'  (psi)’.3X. 

‘Standard  Dsviation  -'.F12.4,'  (psi)’) 

(//IX. ‘This  is  ths  END  I ’//) 

FORMAT  (/IX.'Movs  taskmastsr  again  and  toks  doto?  (y/n)  [DEF  -  Y]'/) 
FORMAT  ( 1X.2F1 1 .4) 

FORMAT  (1X.2F7.2.F102.F8.2.F10.2.F8.2.1P3E12.4) 

FORMAT  (II) 

FORMAT (1X.2F7.2.F10.2.8X.F10.2.8X.1 P3E 12.4) 

F0RMAT(1X. 13, 1P10E12 . 4) 

FORMAT (IX. ‘CRA. DAT.  MACH  NO  -  '.F4.2,'.  FREE  STREAM  VELOCITY  -  ‘.F5.1. 
k  ‘.  ALPHA-  ‘ .F4. 1 ./IX. 'Ssptsmbsr  9.  1987' .4*. i4. ‘  hrs.‘ 

k  //4X.  MROF  NSAkFL  NCOEF'/7X.  I2.6X.  I2.6X.  I2/1X.  ‘  ID 

k  Z-POS  U  V  W 

*  Vt2  Wt2  UV  UW‘ ) 


X-POS', 
Ut2‘ . 
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SUBROUTINE  S£TUP(CLK_lO,OEV_lO.RATE.MCHANL.t«*T,l*T$) 


Subroutine  written  by  Moj  Richard  R  Ryles,  NPS  Student,  13  Moy  87 

Updoted  15  Jul  87 

This  subroutine  is  intended  to  attach  the  internal  clock  ond  A-0  board 
to  the  wain  program  and  allow  them  to  be  used  for  data  acquisition. 

This  subroutine  should  be  used  with  subroutines  GETDATA.  DETACH 


INCLUOE  SYSSLIBRARYLIOSET.FOR’ 


Declare  local  variables 


CHARACTER* 4  CLK.DEV 
CHARACTER* 1  ANS 

INTEGER  SYS.STAT 

IStatus  ret 

by  LIO  calls 

INTEGER  0EV_I0 ,CLK_ID 

ILIO  device 

ID 

INTEGER  NCHAN 

IChannel  f 

INTEGER  NPT. NPTS 

REAL  RATE. TIME 

Input  desired  clock  and  somple  parameters 


WRI TE( • . 1 ) 

REAO(* .2)NCHANL 
RRITE( • .3) 

READ( • . • )RATE 
IF(RATE  GT . 600 . 80) THEN 
*IRITE(*  .4) 

GOTO  10 

ELSE 

IF (RATE. LT.0.01)WRIT£(*,5) 
IF(RATE.LT.0  01)GOT0  10 

ENOIF 

RRITE(* .6) 

READ(* .7)  NPT 
I F(NPT . GT . 5000) THEN 
#RITE(*.8) 

GOTO  20 

ELSE 

TF(NPT.LE  0)WRITE(*,9) 

I F(NPT  LE . 0)GOTO  20 

ENOIF 

NPT  ■  NPT  ♦  1 

TIME  -  FLOAT (NPT- 1 )/RATE 

»IRITE(  •  .  1 1 )  TIME,  TIME/60. 

NPTS  -  2  •  NPT 
CLK_DEV  -  ' KZA0 ' 

*RITE(*. 12) 

REAO( • . 1 3)ANS 

IF  (ANS.EQ  'y-  OR  ANS.EO. -Y’)  THEN 
WRITE(»,14) 

READ(* , 15)CLK_DEV 


Attach  the  KWV11-C  Clock  and  set  up  for  mapped  I/O 


SYS  STAT  -  LIOSATTACH  (CLK  ID,  CLK.DCV.  LlO$K_MAP) 

IF(.NOT  SYS  STAT)  CALL  L1BSSIGNAL(XVAL(SYS.STAT)) 
SYS.STAT  -  L IOSSET_R  (CLK  10,  LIO$K_CLK_RATE ,  1 .RATE) 

I F (  NOT  SYS  STAT)  CALL  LIBSSIGNAL(*VAL(SYS_STAT)) 
SYS  STAT  -  LIOSSET  I  (CLK  ID.  LIO$K_TRIG.  1.  LIO$K_II*IEDIATE) 
I F (  NOT  SYS.STAT)  CALL  L IB$SIGNAL(XVAL(SYS_STAT ) ) 


Attach  the  ADVil-0  ond  set  up  for  mopped  I/O 


SYS.STAT  -  LIOSATTACH  (0EV_ID,  'A2A0-,  LIOSK.MAP) 


pv.www.iw  "A-v.  7rx*?r  v*  ’>  •>'  ■ ,» vv  *v  "v  >  \j*  vv  i^v’vwuv^vvi.'vi.^'.’vw'.'vim.v.v  v*  wn.'iL'iu'v.'i'  k-  .<_-  <- 


ir(  .MOT.  SYS.STAT)  CALL  LlB$SlGHAL(XVAL(SYS_STAT)) 


C  S«t  up  the  AOV  for  synchronous  tronsfer 
C  Sot  channel  (nehonl) 

C  Gain  of  1 

C  Trigger  on  LIOSREAD  and  fill  buffor  at  fast  at  pottibio 

C - — - - - - 

SYS.STAT  -  LIOSSET.I  (0EV.10,  L!0$K_AD_CHAN .  1 , NCHAML . ) 

IF(.NOT.  SYS.STAT)  CALL  LI8SSIGNAL(XVAL(SYS  STAT)) 

SYS.STAT  -  LIOSSET.I  (0EV.1D,  L I OSK.AD.GA I N ,  1.1.) 

IF(.NOT.  SYS.STAT)  CALL  LIB$SIGNAL(XVAL(SYS  STAT)) 

SYS.STAT  -  L I OSSET. I  (OEV  10.  LIOSK.TR !G . 2 . LIOSK  CLK  POINT, CLK  ID) 

1F(  NOT  SYS.STAT)  CALL  LIB$S1GNAL(XVAL(SYS  STAT)) 

c - 

C  Format  Statements 

C - 

1  F0RMAT( 1 x ,' Enter  chonnol  number  to  sample:  (8, 1.8. 9)’) 

2  FORMAT (12) 

3  FORMAT (/I*,’ In put  desired  Clock  Rato  ( 608  Hz  -  0.01  Hz):'/) 

♦  FORMAT ( 1 x , ' Fo I  I ow i ng  directions!  Input  too  h i gh . ’//20x . ' Try  again  '/) 

5  FORMAT ( 1 > , ' Too  low!  Lots  of  time  on  your  hands? ‘//20x . • T ry  again  '/) 

6  FORMAT(/1 x . ’ No .  of  data  points:  (any  number  less  than  5,001)’/) 

7  FORMAT (14) 

8  FORMAT (/I x Hey !  Follow  directions!  Too  high.  Try  again!  '/) 

9  FORMAT (/I x Hey !  Last  input  is  nonsense!  Try  again!’/) 

11  F0RMAT(/1 x . ’ The  approximate  run  time  for  eoch  port  will  be  ’  .  f  S  .  1 . 

k  ’  seconds  (or  '.f5.2.‘  minutes)'/) 

12  FORMAT(/1 x . ' Do  you  wont  to  change  the  Clock  Device  Nome  from  [KZA0]?’, 

*  •  (y/n)  (DEF  -  no]’/) 

13  FORMAT ( A1 ) 

14  FORMAT(1x, 'Enter  Clock  Device  Nome:  (K2A0 , KZB0) ' ) 

15  FORMAT (A4) 

C - - - — 

RETURN 

EfC 
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SUBROUTINE  GETOATA(DEV_IO.f#»T.NPTS, VOLTAGES) 

Subroutine  written  by  Mo j  Richard  R  Rylee,  W>5  Student,  28  Apr  87 

Update  IS  Jul  87 

This  subroutine  it  intended  to  acquire  the  data  from  the  A-0  board 
with  the  specified  parameters  from  subroutine  SETUP.  It  also  will 
translate  the  A-0  data  into  voltages  to  be  converted  to  physical 
quantities  in  the  main  program. 

This  subroutine  should  be  used  with  subroutines  SETUP .DETACH . OUTPUTD 
INCLUDE  ' SYSSLIBRARY : L 10SET . FOR ' 

Declare  local  variables 

INTEGER  SYS.STAT  iStatus  ret  by  LIO  calls 

INTEGER  OEV.IO  ! LIO  device  10 

INTEGER  DATA_LENGTH  INumber  of  bytes  of  data  read 

INTEGER  NPT.NPTS  ! Number  of  points  to  sample 

Declare  dota  buffers 

INTEGERS  RAW_OATA(5000)  15000  point  raw  data  buffer 

REAL  VOLTAGES (5000)  15000  voltages  output 

Get  a  raw_dota  buffer  of  'npt’  values 
Uses  LIOSREAD  to  reod  the  A/D 
Note  length  of  buffer  is  in  bytes 
As  Is  the  returned  data  length 

SYS  STAT  -  LIO$READ(DEV  ID.  RAW  DATA,  NPTS,  DATA  LENGTH.  ) 

I F( . NOT .  SYS_STAT )  CALL  LIB$S1GNAL(XVAL(SYS_STAT) ) 

Convert  the  raw  data  to  voltages  using  LSP$FORMAT_tranSLATE_ADC 

CALL  LSP$F0RMAT.TRANSLATE_ADC(RAW_DATA,  VOLTAGES,  NPT ,  ,  .  ) 

RETURN 

END 


o  o  o  o  o  o  o 


subroutine  OETACH(CLK_IO.OEV_10) 


Subroutine  written  by  ktoj  Richard  R  Ryles.  NPS  Student  26  Apr  87 

Update  16  jun  67 

This  subroutine  is  intended  to  detach  the  internal  clock  and  A— 0  board 
from  the  dam  program  and  allow  them  to  be  used  for  others  application 

C  This  subroutine  should  be  used  with  subroutine  SETUP 


INCLUDE  'SYSSLIBRARY  LIOSET  FOR' 

c*  • » •  • . . . . . . 

C  Declare  local  variables 

C* . . . . . . . . 

INTEGER  SYS.STAT  ‘Status  ret  by  LIO  calls 

INTEGER  DEV  IO.CLK  ID  *  L 10  device  ID 

. . . . 

C  Oetach  from  the  Clock 

. . * . * . . . . . 

SYS_STAT  -  LIOSDETACH  (CLK_ID  ,  ) 

IF(  not  SYS  STAT)  CALL  LIBSS!GNAL(XVAL(SYS_STAT)) 


SUBROUTINE  OUTPUTD(RATE. NPT. NCHANL.  AMEAN,  $TD_DEV.J.VOt.TACES) 


Subroutine  written  by  Uoj  Richard  R  Rylee,  NPS  Student  28  Apr  87 

Update  24  Jui  87 

This  subroutine  is  intended  to  output  diagnostic  results  to  o  file. 


This  subroutine  should  be  used  with  subroutine  SETUP,  DETACH.  GETDATA 


Dec  lore 

locol  variobles 

1 

INTEGER 

NCHANL 

(Channel  f 

INTEGER 

J 

(General  loop  indexes 

INTEGER 

NPT 

(Number  of  points 

REAL*8 

AMEAN, STD  DEV 

REAL 

RATE.VOLTAC£S(5000) 

Output  dato  to  a  file 


OPEN (UN IT  -  10,  FILE  -  1 SCANIDI AG  DAT ‘ .  STATUS  -  'NER' ) 

WRITE( 10,1)  J.  RATE,  NPT-1 ,  NCHANI ,  AMEAN ,  STD.DEV 

FORMAT(/1 x . ’ Scon i va I v«  § ' , i 3/1 x , 'C I ock  Rote  , f 6 . 1 ,5* , ’ Do  to  Points', 
‘Sampled  . i 4/1 *, 'Channel 5* . 'Mean ’,4*. 'Standard  Deviation'/ 
3*. i 3 . 2f 12.4/) 

NR  I TE( 10 . • ) ' Scan i vo I ve  f’.J 
DO  20  1-1 , NPT 

NRITE(10,») ' I-' . I ,VOLTAGES(I ) 

CLOSE( 10) 

RETURN 

END 


on  ooooo  o  o  o  o  o  ooooooooooooooooooooooooooo  o  o 


SU8A0UTINE  *W*A(NSMf>LES,  ANSWER.  8LUE_GREEN) 
PURPOSE 

Drives  the  ocquieition  o f  the  LDV  data 
ERROR  HANOI INC 

IER  •  1  meons  no  'fatal"  error. 


EXTERNAL  references 

NAME  DESCRIPTION  AND  SOURCE 

DVWA  Calls  DRV11-WA  dsvics  drive r  to  acquire  A/D  dota 


DEVELOPMENT  HISTORY 


DATE 

INITI, 

07/26/83 

TML 

12/12/85 

CLH 

02/25/87 

GBG 

05/06/87 

CLH 

06/08/87 

RRR 

07/15/87 

RRR 

07/24/87 

RRR 

08/29/87 

RRR 

DESCRIPTION 
Adapted  from  OAPNT 
Modified  to  print  results  in  either 
octal  or  i n tege r 

Modified  for  use  with  0RV11-WA  and  VMS 
Modified  to  accept  "nee"  DAMUX1 . 

Modified  to  subroutine  *  deleted  CIH  S/87  Mod 

Modified  to  calculate  velocities  +  optimized 

for  speed  of  execution 

Modified  to  vie*  Histogram 

Modified  to  occept  calibration  constant 


AUThOR(S):  Ted  Lichtenstein,  Informatics  General  Corp. 


IMPLICIT  NONE 


PARAMETER  MtAWMX-80000 

PARAMETER  IUNTI-5 

PARAMETER  LUNOUT-6 


INTEGERS 

INTEGERS 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

REAL 

REAL 

REAL 

LOGICAL 

CHARACTER*  19 

CHARACTER*  1 


I RAWDT ( NRAWMX ) 

1 1 RAWDT ( NRAWMX ) , 1 2RAW0T (NRAWMX ) 

IER.HlSTO.Nf r. ONE, TWO 
IYESNO,  NRAWIN.  I,  BLUE.CREEN 
NCHANS ,  NSAMPLES,  INOEX,  JNDEX 
NS . ANSWER . COUNT , END 

TIME(NRAWMX) , FREQ(NRAWMX) . VELOC ( NRAWMX ) . I 3RAWDT ( NRANMX ) 

MIX(2) , GAM(2) ,DEL(2),KAP(2) .D(2) .1(2) ,MEAN( 2 ) . FREOUENCY ( 2 . 8001  ) 
FREQONE . FREQTWO , FRE(NRA*MX ) . FREQ1 _AVE . FREQ2  AVE . CAL_C0NST ( 2 ) 
DATSELCT  (3) 

DIRECTION 

ANS 


DATA  MIX(1 ) . M I  X ( 2 )  /  -1.E-8.-1.E-6/ 
INCLUDE  ' (SSYSSRVNAM) ' 

100  CONTINUE 


Fill  the  raw  data  buffer  with  this  bit  pattern: 
"1010101010101010" 

This  equals  125252(8).  -21846(10). 


DO  200  I  -  1 ,  NRAWMX 

I RAWDT ( I )  -  '  1 25252 ' 0 
hrawot(i)  -  '000000'o 

1 2RAWDT ( I )  -  ' 000000  o 
! 3RAWDT ( I )  -00 
200  CONTINUE 


C 


Input  number  of  samples  to  be  ocquired 
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NCHANS  -  4 

300  WRITE  (iunti,  ie) 

10  FORMAT  ('  Enter  number  of  somples:*) 

READ  (LUNTI,  .)  NSAAffLES 
NR  AW  IN  -  NSAAPLES  •  NCHANS 

IF  (NRAWIN  GT.  32767  OR.  NRAWIN  . LE .  0)  THEN 
WRITE  (LUNTI .  20) 

20  FORMAT  ('  The  number  of  samples  •  number  of  channels  must  be  32.767' 

k  *  or  I  ess . ' ) 

GO  TO  300 
ENO  IF 


Set  up  the  LDV  parameters  prior  to  acquisition 


I F(8LUE_GREEN . EQ . 1 ) THEN 
GAM( 1 )  -  0  488 
GAM( 2)  -  0  5145 
IF(CAl_CONST(1).EQ.1.)  THEN 

WRI TE( *.«)’ Enter  Blue  calibration  constont’ 

READ( • . • )  CAL_C0NST(1) 

END  IF 

I F (CAL_C0NST ( 2 )  EO . 1 . )  THEN 

WRITE( Enter  Green  calibration  constant' 

READ( • . • )  CAL_C0NST(2) 

END  IF 
ELSE 

GAM(1)  -  0.5145 

GAM( 2 )  «  0.488 

I F (CAL_CONST ( 1 ) . EQ. 1 . )  THEN 

WRITE( •,•)' Enter  Green  calibration  constant' 

READ( •  .  • )  CAL_CONST ( 1 ) 

END  IF 

I F ( CAL_CONST ( 2 ) . EQ . 1 . )  THEN 

WRITE(» ■•)' Enter  Blue  calibration  constant' 

READ( • . • )  CAL_CONST (2 ) 

END  IF 
END  IF 

IF(MIX( 1 ) . LT . 0.0)  THEN 

I F(BLUE_GREEN  EQ . 1 )WR I TE( • . • ) ' Input  BLUE  mixing  frequency  (MHz)' 

I F(BLUE_GREEN . EQ. 0)WRITE( • . • ) ' Input  GREEN  mixing  frequency  (MHz)' 
READ( • . •)  MI X( 1 ) 

END  IF 

I F(MI X( 2)  LT  0  0)  THEN 

I F(BLUE_GREEN . EQ . 1 )WRITE(» . • ) ' Input  GREEN  mixing  frequency  (MHz)’ 
I F(8LUE_GREEN . EQ.0)WRITE( • . •) ' Input  BLUE  mixing  frequency  (MHz)’ 
READ( • , • )  MI X( 2 ) 

END  IF 

I F( ANSWER  EO  0)  THEN 

WR I TE (•■•)' Do  you  wont  probe  volume  specifications?  (y/n) ' 

REAO( • . '(a)')  ANS 

I F ( ANS  EO  ' Y '  OR.ANS.EQ. 'y')  ANSWER  -  I 
END  IF 

WR I T  E ( • . • ) '  ' 


Select  the  data  to  be  returned. 


DATSELCT  (1)  -  FALSE. 
DATSELCT  (2)  -  FALSE. 
DATSELCT  (3)  -  TRUE. 


LDV  Doto  Acquisition 


ONE  -  1 
TWO  -  2 

0PEN(UNIT-5,FILE-'VEL0CITY.DAT' ,STATUS-'NEW' ) 

DO  INDEX  -  1.  NSAMPLES 

CALL  DAMUX  (NCHANS.  ONE,  DATSELCT.  .  .  IRAWOT,  IER) 
!F(IER  NE  I)  WR I TE( • . • ) '  IER  in  DAMUX  -'.IER 
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Conversion  to  Frequency 


DO  1-1.2 

IRAWOT(I)  .  INOT(IRAWOT(I)) 

I3RAW0T(I)  -  0.0 
I 2RAWDT ( 1 )  -  0 
I 1RAW0T ( I )  -  0 

IIRAWOT(I)  -  IIBITS(IRAWOT(1).10,4) 

I 2RAWDT ( I )  -  I IBI TS( IRAWDT ( I ) , 0 , 1 0) 

I 3RAWDT ( I )  -  (2 . I IRAWDT ( I ) )  •(  f loot ( I2RAWDT( I )) ) 
TIME(l)  -  I 3RAWDT ( I )/32 . 

FR£0( I )  -  1000 . /TIM£( I ) 

END  DO 

FREQUENCY ( 1 , INDEX)  -  FREQ( 1 ) 

FREQUENCY (2,1 NOEX )  -  FREQ(2) 


Convert  frequency  dato  to  velocities 


CALL  VELOCITY  (FREQ.  VELOC  ,  TWO  ,  MI  X  .ANSWER  .  GAM.  DEL  ,KAP .  D ,  I  .  Nf  r  .  CAL.CONST  ) 
WRITE(5.30)  VELOC ( 1 ) , VEL0C(2 ) 

FORMAT ( 1 X , 2F1 1 .4) 

END  DO 
CLOSE(5) 


Calculate  average  counter  frequency  and  print  to  sceen  for  diagnostics 


FREQONE  -  0  0 
FREQTWO  -  0.0 
DO  34  I  -  I.NSAfcPLES 

FREQONE  -  FREQONE  +  FREQU£NCY( 1,1) 

FREQTWO  -  FREQTWO  +  FREQUENCY (2 , I ) 

FRE01  AVE  -  FREQONE/FLOAT(NSAkPLES) 

FRE02_AVE  -  FREQTWO/FLOAT(NSA»*>LES) 
nKRiTE(..«r  • 

WRITE(».«)  Average  counter  frequency  f  1  -',FREQ1_AVE 
WfR  I  TE(«  .•)  'Average  counter  frequency  #  2  »',FREQ2_AVE 


Write  Probe  Volume  dimensiont  and  parameters  to  a  file 


I F( ANSWER  EQ  1 ) THEN 

OPEN(UNIT-1 5. FILE- •Probe. dot ' . STATUS- ' NEW ' ) 

WRITE( 1 5 , • ) 'Probe  Volume  Data’ 

WRI TE(  1  5 .  •  )  ' - ’ 

k  • - • 

DO  40  I  -  1,2 

I F(BLUE_GREEN . EQ . 1 )  DIRECTION  -  ’Horizontal  Velocity' 

IF(BLUE_GREEN. EQ.0)  DIRECTION-  'Vertical  Velocity 
I F( I . EQ . 2 )  THEN 

I F(BLUE_GREEN . EQ. 1 )  DIRECTION  -  'Vertical  Velocity 
IF(BLUE  GREEN. EQ.0)  DIRECTION  -  'Horizontal  Velocity' 

END  I F 

WRITEC 15,1 )GAM(!) . DIRECT  I  ON , DEL( I ) ,KAP( I ) ,D( I ) , I (I)*  001 .Nf  r 
FORMAT(/1 X . 'Wove  I ength  of  light  <microns>  -  '.f7.4,'  (Measuring', 
k  at 9 , ' ) ’//I x , ' Fr i nge  Spacing  <microns>  -  ' . f 1 0 . 2//1 x , ' Lens '  . 

k  '  angle  <deg>  •' . f 7 . 2//1 x . 'Probe  volume  diameter  <microns>  - 

k  f 10. 5//1 x Length  of  measuring  volume  <mm>  -  '.f10.4//1x, 

k  Number  of  Fringes  -' . i 3/ ) 

WRITE( 15. •) ' - 

k  • - • 

CONTINUE 
ClOSE( 1 5 ) 

END  IF 


Option  to  print  o  histogram  of  velocity  doto  on  screen  for  diagnostics 


WR I TE( • , • ) '  ' 

WR I TE( • , • ) ' Do  you  -ant  o  Histogram?  (1/0)' 


o  o  o 


READ(» . •)  HISTO 
CNO  —  NSAtflES 

If(HISTO.EQ. 1 )CALL  HISTOCRAM(COUNT  . END . BLUE.GREEN) 


Ago i n? 


WRITE( LUNTI , •) '  1 

WRITE  (LUNTI.*)  'Acquire  LDV  Data  again'’  (0-NO,  1-YES)  ’ 
READ  (LUNTI,.)  IYESNO 
WR I TE( »  ,  » ) '  1 

IE  (  IYESNO. EQ. 1  )  GO  TO  100 

RETURN 

END 
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SUBROUTINE  VELOCITY(FREQ,VELOC.N. MIX.  ANSWER. GAM. DEL. KAP.  to.  I  ,Nf  r. 

k  CALJDONST) 


REAL  THETA.  GAMMA,  DELTA,  d,  f  .  E  .  KAPPA ,  Dsubm .  I  m .  De .  FREQ(N)  .  VELOC(N)  ,GAM(  2) 
REAL  SHIFT. MIX(2).0EL(2). CAL  CONST (2) . KAP(2) . to(2) . 1(2) 

INTEGER  I ANS. ANSWER,  Nfr 

if(cal_const(i)  eo  0  0)  cal_const(i)  -  i  0 
I F(CAL_C0NST(2)  EQ. 0.0)  CAL_C0NST(2)  -  10 
PI  -  ASIN( 1 .0) »2 . 


d  -  13.  ’  Bern  spacing  (im) 

f  •  600.  I  Focal  length  of  transmitting  lens  (im 

E  ■  1.  I  Beam  expansion  ratio 

KAPPA  -  ATA*C(E*d/2 . /f )  I  Lens  angle  (deg) 

00  5  I  -  1  .2 

GAMMA  m  GAM( I )  !  Wavelength  of  Laser  (microns) 

OELTA  -  CAL_C0NST ( I )  .  GAMMA  /  (2 . »SIND(KAPPA) )  !  Fringe  spacing 
De  »  1  .  1  !  Diameter  of  the  incident  laser  beam  (im 

Osubm  «  (4 . *GAMMA« f/(PI «De«E) )  !  Probe  voluse  diometer  (microns) 

Im  »  Dsubm  /  TANO(KAPPA)  !  Probe  volume  length  (microns) 

Nfr  ■  1  27  •  d  /  De  1  Number  of  fringes 


DEL(I)  -  OELTA  !  V 

KAP(I)  -  KAPPA  ! 

Dm( I )  »  Osubm  ! 

1(1)  -  Im  ! 

SHIFT  -  40.0 

FREO(I)  -  FREQ(!)  -  SHIFT  ♦  MIX(I) 
VELOC(i)  -  OELTA  .  FREO(i) 

RETURN 

ENO 


Probe  volume  length  (microns) 

Number  of  fringes 

Variables  possed  to  output  file 


•  Bragg  shift  (MHz) 

MIX  «  Mixing  Freq  (MHz) 

FREQ  -  Adjusted  Doppler  Freq 

(MHO 


risVJ 


•Vi 


V  v"  **  w* V*  v*  'V 


>>W/V 


o  o  o  ooo*  o;  non  w  —  o  o  o  n  nnnnnnnn 


SUBROUTINE  H I  STOGRAM(  COUNT  ER.E)©,BCUC_  GREEN) 

Subroutine  or  it  ton  by  Mo j  Richard  R  Ryles.  M>S  Student,  it  Jun  87 

Update  29  Aug  87 

Thie  eubroutine  it  intended  to  creote  too  histogram*  of  velocity 
dato  from  a  too  component  LDV. 


CHARACTER *4  STRING 

CHARACTER*22  LABE LX  I  variable*  for  the  exit 

CHARACTER* 17  LA8ELY  I  lob  lee 

CHARACTER  *29  LAB  ELM 

INTEGER  STATUS, I. I SHADE. COUNTER. EM) 

INTECER  N(3) ,BLUE_GREEN 

REAL  XCONTROLI4)  .XMIN,  XMAX 

REAL  YC0NTR0L(4),CCUNT(8) 

REAL  XLCW(8) 

REAL  XHIGH(S) 

REAL  Y(8) .VELOCITY (8288, 2) ,X(828t) 


Nt 

>1) 

-  i 

N( 

2 

-  i 

N( 

0  0  9 

!s) 

• 

-  i 

Rood  Velocity  Data  Tile 


I  number  of  doto  point* 

I  number  of  column*  of  data 
I  number  of  roo*  of  data 


OP EN(UNIT-5. FILE-’ VELOCITY. DAT’  , STATUS- ’ OLD ’  ) 
IF(BLUE  GREEN. EQ. 1) THEN 
DO  18  I  -  1 .END 

READ(5 , t )  VELOClTY(!.1).VELOCITY(I ,2) 

ELSE 

DO  28  1-1 .ENO 

READ(5 . 1 )  VELOCITY( 1,2) , VEL0CITY( 1.1) 

END  IF 


Find  the  Win  and  Max  values 


CL0SE(5) 

DO  88  J  -  1  .2 

DO  38  K  -  1 .end 

X(K)  -  VELOCITY  (K.j) 
XMIN  -  1888.8 
XMAX  -  -1888.8 
DO  48  I  «  1 .end 

IFfXMIN. GT . X( I ) )  XMIN  -  X(I) 
I F(XMAX . LT . X( I ) )  XMAX  -  X(I) 

Initialize  X-Axi*  controls 


XCONTROLf 1 )  -  5.5  !  length  of  x-oxi*.  in  inches 

XCONTROC(2)  -  XMIN  !  left  most  value  on  x-axis 

XCONTROLl 3)  -  XMAX  !  right  most  value  on  x-axis 

XCONTROL(4)  -  (XMAX-XM!N)/8.  I  increment  of  plotted  points 


Find  y  dato  from  sample  array 


DO  58  I  -  1.8 
Y(  I )  -  8.8 

XLOdf(I)  -  XMIN  +  FLOAT (1-1  )«XCONTROL(4) 
IF(I.EQ.8)THEN 
XHIGH(S)  -  XMAX 

ELSE 

XHICH(I)-  XMIN  ♦  FLOAT(I)*XCONTROL(4) 

END  IF 

58  COUNT ( I )  -  XMIN  ♦  FL0AT(I)*XC0NTR0L(4) 

C 

XTOTAL  -  88 


o  o  r>  o  o  o 


60 

C 


70 


DO  00  1  -  1  .and 

XTOTAL  -  X(I)  ♦  XTOTAL 

lF(XM)  .GE.XMIN.AM).X(I)  . LT .COUNT( 1 ) )  Y(1)  -  Y(1)  ♦  1.0 
IFfXi I ] . GE. COUNT (1 ) .ANO .X( 1 ) . IT. COUNT (2))  Y(2)  -  Y(2)  >  1 
IFtXf I l GE. COUNT! 2) ANO.Xi  I ) . LT. COUNT (3))  Y(3)  -  Y(3)  ♦  1 
IFfXf 1 1  GE .COUNT! 3) . ANO. XI I  ) . LT .COUNT l4l)  Y  4  -  Y<4)  ♦  1 
IFi Xl I ) .GE .COUNT I  *) . AND . XI I ) . LT .COUNT (3) I  Y(5)  -  Y(5)  1 

IFiXi 1 1  GE. COUNT (5) . AND . Xl  I ) . LT .COUNT (0)1  y(0)  >  Y(0)  ♦  1 
I F I X ( 1 1 . GE . COLNT  1 0  > .  AND .  Xt  1 1  .  LT  . COUNT (  7  )  )  Y(7)  -  Y(7)  ♦  1 
1F( Xl  I ) GE. COUNT (7) .ANO.XI I ) . LT .COUNT(8) )  Y(8)  -  Y(8)  ♦  1 
IF(X(1)  GE.COUNT(8)  ANO.X(I)  .IT.XMAX)  Y(8)  -  Y(8)  ♦  1.8 
XMEAN  -  XTOTAL/FLGAT(ENO) 

YUAX  -  0.0 

DO  70  I  -  1.8 

Y(l)  -  Y(  I )/FLOAT (DC) 

I  F(YI4AX .  LE . Y( I )  )  YMAX  -  Y(I) 

Y( I )  -  100.  •  Y(I) 

«IR1TE(*,2)YUAX»100.  .XMEAN 


Initial i z •  Y  axis  control* 


YCONTROL(I)  -  4.5 
YCONTROL(2 )  -  0.0 
YCONTROL(3)  -  YMAX*  1 20 . 
YCONTROL(4)  -  5.0 


Histogram  lab«l* 


I  height  of  graph  in  inch** 

I  lo«*r  boundary  valu* 

I  upp*r  boundary  valu* 

I  incr*i**nt  of  plotted  point* 


LABELX  »  '  VELOCITIES  (mp*)  !  x— axis  lob* I 

l F(COUNTER.GT.0) LABELX  -  ‘FINAL  VELOCITIES  («p*)‘  I  x-axi* 

LABELY  -  PERCENTAGE  OF  PTS‘  I  y-axi*  lab* I 

LABELM  -  'HORIZONTAL  VELOCITY  HISTOCRAM'  !  uppar  lab* I 

IF(J.EQ.2)  LABELM  -  ‘  VERTICAL  VELOCITY  HISTOGRAM- 
STRING  -  ’ IXSY '  I  no d*  string 

I SHADE  >2  !  filling  of  graph 


C  S«t  up  coordinate  *y*t*m  with  valu** 

. . . . 

IF(XMIN.GE.XMAX)GOTO  90 

CALL  LGP$PL0T(1 .STRING. .0.N. LABCLX. LABELY. STATUS. . . 

X CONTROL , Y CONTROL , . LABELM) 

IF  (STATUS  .NE  LGPS_SUCCESS)  THEN 
CALL  LIB$SIGNAL(XVAL(STATUS) ) 

ENDIF 

. . . 

C  Call  rout  in*  and  plot  histogram 

C  *  *  •  •  • . . 

CALL  LGP$HIST(1 .XLOW.XHIGH.Y. 8. STATUS, .ISMADE) 

IF  (STATUS  NE  LGP$_SUCCESS)  THEN 
CALL  LIB$SIGNAL(XVAL(STATUS)) 

ENDIF 

PAUSE  'Raiting  for  v i «vi ng/pr int i ng  Histogram;  Typ*  CONT  »h *n  finished' 
CALL  LGP$TERMINATE_PLOT(l)  I  t*n*inat*  plotting 

80  CONTINUE 
return 

90  RRITE(*.*)'  ' 

»FITE(*. •) '•••  No  Histogram  created;  xmin>«xhigh  •••' 

GOTO  80 


C  Formats 

. . ♦  •  •  • 

1  FORMAT ( 1X.2F1 1  .  ♦) 

2  FORMAT (/IX. 'MAXIMUM  PERCENTAGE  OF  SECTORS  -  '  .  F6  1//1X .  '  SAAFLE  MEAN  -  '. 

0  F9.2) 

run 


lab*  I 


P 


SUBROUTINE  REDUCE_DATA(NSAMP. COUNTER. ENO. MEAN. STD.DEV1 .STD  DEV2. 
k  BLUE.GREEN) 

C - 

C  Subroutine  written  by  Moj  Richard  R  Ryle*,  NPS  Student,  12  Aug  67 
C  Update  29  Aug  87 

C 

C  This  subroutine  is  excludes  all  velocity  data  from  the  LOV  thot  fall 

C  outside  two  standard  deviations  from  the  mean. 

c - 

REAL*8  ST0_0EV1 , STD_DEV2 ,NEW_STD(2 ).TOTAL(2). SUM (2 )  ,  X  ,  Y 
REAL  VELOC (2.1 0000 ) . MEAN ( 2 ) , UP , DOWN , TOLERANC E 
REAL  NEW(2)  ,VELO(2, 20000) 

INTEGER  NSAMP .COUNT (2) , TOT , END , BLUE_GREEN 

. . .  •  . . . 

C  Read  velocity  data  from  file 

C..  . . . 

OPEN (UN  IT-3,  FILE- 'VELOCITY. DAT1 , STATUS- ' OLD ' ) 

00  1 0  I  -  1  .  NSAMP 

10  REA0(3 , 1 ) VELOC ( 1,1) .VELOC (2 . I ) 

CL0SE(3) 

. . . . . 

C  Calculate  range  of  acceptable  values  i.e.  all  velocity 
C  data  that  falls  within  two  standard  deviations  of  the  mean 

. . . . 

DO  30  J  -  1 . 2 

IF( J . EO. 1 )THEN 

TOLERANCE  -  ST0_DEV1 *2 . 

UP  -  MEAN( 1 )  +  TOLERANCE 
DOWN  -  MEAN( 1 )  -  TOLERANCE 

ENO  IF 

I F( J . EO . 2) THEN 

TOLERANCE  -  STD_DEV2»2. 

UP  -  MEAN(2)  ♦  TOLERANCE 
DOWN  •  MEAN ( 2 )  -  TOLERANCE 

ENDIF 

. . . . 

C  Flag  each  point  outside  of  two  standard  deviations 

C  •  •  •  • . . 

DO  20  I  -  1  .NSAMP 

20  I F (VELOC( J . I ) . LT . DOWN . OR . VELOC ( J , I ) .GT . UP)VELOC( J , I )— 1000  0 

30  CONTINUE 

C.  . . . . . 

C  Calculate  new  mean  and  standard  deviations 

. . . . . 

DO  60  J  -  1 . 2 

TOTAL(J)  -  0.D0 
COUNT ( J )  -  0 
SUM(J)  -  0  D0 
DO  50  I  -  1  .NSAfcP 

IF(VEL0C(1 . I)  .NE. 1000.0. AND. VEL0C(2. I)  NE. 1000.0)  THEN 
TOTAL(J)  -  TOTAL(J)  +  VELOC(J.I) 

COUNT ( J )  -  COUNT ( J )  +  1 
ENDIF 

50  CONTINUE 

I F ( COUNT ( J ) . EO -  0 )  COUNT ( J )  -  1 
NEW(J)  -  TOTAL(J)/FLOAT(COUNT(J)) 

DO  55  K  -  1  .NSAMP 

I F (VELOC( 1 ,K)  NE. 1000,0  AND. VELOC (2, K)  NE  1000  0)  THEN 
X  -  VELOC(J.K)  -  NEW(J) 

Y  -  X.X. ( 1  /FLOAT (NSAMP- 1 ) ) 

SUM(J)  -  SUM(J)  +  Y 
ENO  IF 


55  CONTINUE  j 

NEW_STD( J )  -  SORT ( SUM( J ) )  i 

60  CONTINUE  ‘ 

. . . .  ; 

C  Output  new  values  to  screen  .  “ 

C  •  •  . . .  , 


76 

f 

t 


WRITER,.)'  • 

WRITE(» , •) ’New  cleoned-up  data’ 

IF(BLU£_GREEN.EQ. 1)TH£N 

WR I  TE(  •  .  • )  ‘  U  dir’,NEW(1).NEW_STD(1) 

WR  I  TE(  •  ,  » )  ’  V  dir’  .N£W< 2 ) . NEW_STD( 2 ) 

ELSE 

WR I TE( • . • ) 'U  dir’ ,  NEW(2) , N£W_STD(2) 

WR ITE( • . • ) ' V  dir' ,NEW( 1 ) , NEW_STD( 1 ) 

END  I F 

MEAN( 1 )  =  NEW( 1 ) 

MEAN(2)  -  NEW(2) 

STD  DEVI  -  N£W_STD( 1 ) 

STO  DEV2  -  NEW  STD{ 2 ) 

. . . . .  . . . 

C  Write  interim  data  file  containing  just  good  data 

. . . . .  •  . . 

OPEN( UN  IT-5,  FILE-' INTERIM.DAT ’ . STATUS- ’ NEW’ ) 

DO  70  J  -  1  ,2 

DO  70  I  -  1  ,  NS  AMP 

1 F( VELOC( 1 .  I)  N£  1000.0. AND. VELOC(2, I ) . NE . 1000.0)  THEN 
WRITE(5 , 2) VELOC( J . I ) 

END  IF 

70  CONTINUE 

ClOS£(5) 

C  ........  . . . . . 

C  Write  final  velocity  data  file  for  plotting 

C.  . . . . . 

OPEN(UN IT-7, FILE- ’VELOCITY  OAT' , STATUS- ' NEW’ ) 

OPEN( UN  IT-9. FILE-’ INTERIM.DAT’ , STATUS- ’OLD ’ ) 

TOT  -  COUNT ( 1 )  +  COUNT ( 2 ) 

DO  80  I  -  1 .  TOT 

I F( I  LE.COUNT(1))THEN 
READ(9 . 2)  VELO(I.I) 

ELSE 

READ(9 . 2)  VELO(2. I) 

END  IF 
80  CONTINUE 

CLOSE(9 ) 

END  -  COUNT ( 1 ) 

COUNTER  »  1 

DO  1 00  I  -  1  .  END 

100  WR I TE( 7 . 1 )  VELO( 1,1). VELO( 2 , COUNT ( 1 )+I ) 

CLOSE( 7 ) 

RETURN 

C  ........  . . . . . 

C  Format  statements 

. . . . . . . 

1  FORMAT (1X.2F11.4) 

2  FORMAT ( 1 X . F1 1 .4) 

END 
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